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1.  INTRODUCTION 


w 


In  the  design  of  a conformal  slot  array  on  the  surface  of  a 
conducting  cylinder,  the  calculation  of  the  mutual  admittance 
is  a crucial  step,  which  has  been  studied  extensively  in  recent  years. 

In  this  paper,  we  summarize,  in  a handbook  format,  all  of  the  final 
tormulas  of  present  some  typical  numerical  data. 

2.  STATEMENT  OF  PROBLEM 

Referring  to  Figure  1,  two  identical  slots,  circumferential  or 
axial,  are  located  on  the  surface  of  an  infinitely  long  cylinder.  The 
geometrical  parameters  are 

R = radius  of  the  cylinder  (2.1) 

(a,b)  = dimensions  of  the  slot  along  (<t>,z)  directions  (a  is 


the  arc  length  along  the  cylinder)  (2.2) 

(Zq,R(})q)  = center-to-center  distances  between  slots  (2.3) 

"o  -/o  * 

0Q  = tan  ^(Zq/R(|)q)  (2.5) 


The  problem  is  to  determine  the  mutual  admittance  between  these  two 
slots  when  kR  is  large. 

1 


- ^ 4.  ^ 


First  let  us  define  mutual  admittance.  Throughout  this  work 


we  always  assume  that 

ii)  the  slots  are  thin,  and  (2. ha) 

(ii)  their  length  is  roughly  a half-wavelength.  (.l.hh) 

ihen  the  aperture  field  in  each  slot  can  be  adequately  approximated  by 
a simple  cosine  distribution,  which  is  the  so-called  "one-mode”  approximation. 
For  example,  if  slot  1 is  a circumferential  (lower  slot  in  Figure  la), 
its  aperture  field  under  the  "one-mode"  approximation  is  given  by 
(exp  + jwt  time  convention) 


E = , H = I^h^ 


(2.7a) 


where 


1 J ab 
y = R(|i 


(2.7b) 

(2.7c) 


(Vj^,!^)  are  respectively  the  modal  (voltage,  current)  of  slot  1.  The 
mutual  admittance  is  defined  by 


Y = Y 
12  21 


21 


(2.8) 


where  is  the  induced  current  in  slot  2 when  slot  1 is  excited  by 

a voltage  and  slot  2 is  short-circuited.  An  alternative  expression 
for  Y is 


12 


E,  X 14 


(2.9) 


wlierc 


= aperture  of  slot  2 

llj  - in.if’nctic  field  when  slot  1 is  excited  with  voltage  , and 
slot  2 1 covered  by  a perfect  conductor 
= ele(;tric  field  wIhmi  slot  2 is  i;xc  i l ed  witli  voltage  , and 


slot  1 i ; cf)veri:d  by  a perfect  c ondiict  i>;-. 


I 
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•*>  -V  ► 

Because  ” '^9^9’  ^ simple  matter  to  verify  that 

(2.8)  and  (2.9)  are  equivalent  [1]. 

There  is  an  alternative  definition  of  mutual  admittance.  Instead 
of  (2.7),  a modal  voltage  Vj^  (with  a bar)  may  be  defined  through  the 
expression  tor  the  aperture  field  of  slot  1 as  follows: 

E = z V cos  — y (2.10a) 

D I Q 


or  equivalently 
|-b/2 

"l  = 

-b/2 


®)v=0 


(2.10b) 


Then  a different  mutual  admittance  is  defined  by 

(Vj^.V^)  by  (Vj^.V^).  It  can  be  easily  shown  that 

Y = — Y • 

12  2b  12 


(2.9)  after  replacing 


(2.11) 


Two  remarks  are  in  order:  (i)  In  the  limiting  case  that  b -+  0,  Y^.,  goes 

to  zero  as  b , whereas  Y^^  approaches  a constant  independent  of  b. 

(ii)  For  the  special  case  a = A/2  and  R it  is  Y^^  that  is 

identical  to  the  mutual  impedance  between  two  corresponding  dipoles 

calculated  by  the  classical  Carter's  method  [2],  [3],  [4].  (iii)  When 
the  slots  are  excited  by  waveguides  (transmission  lines),  one  often  uses 
Yf2  (Y^.,).  From  here  on,  we  will  concentrate  on  Y^^  instead  of 

The  mutual  admittance  defined  in  (2.8)  and  (2.9)  includes  the  self 
admittance  Y^^^  as  a special  case  which  occurs  when  two  slots  coincide* 
(All  the  formulas  of  Y^^  given  in  this  paper,  except  for  the  one  in 
Section  4,  can  be  used  for  calculating  Y^^  by  setting  i|iy  0 and  z^  -<■  0.) 


i 


4 


3.  EXACT  HUGHES  (GSP)  MODAL  SOLUTION 


Once  the  one-mode  approxim.it  ion  in  (2.7)  is  accepted,  can  be 

determined  exactly  in  te.-ms  of  cylindrical  modal  functions,  as  has  been 
done  by  Stewart,  Golden,  and  Pr idmur^-Brown  [5],  [6].  Tiie  final  result 


reads : 


Circumferential  slots 


“ -.iCnW  z ) 

S ).  )C(m,k,  )e  ^ 

m — '' 


^ ^ .-,1,  (k^b/2)  Alin  (m)i_^  + tt/2)  sin 


2 2 9 

8t!  R (k^l./2) 


(mif  - ir/2) 


(m^j  + tt/2) 

H 


(m^.  - ii/2) 

n 


<>^  = (u/2R) 

cl 


(k  R)  Imk  'r  k R 

C(m,k  ) = Y_  — — k X.  I ? 1 c m t 

0 k.  ..(2),.  . ^ Tur  


•t  ^k^R)  ^k^R/jX 


, if  k 


j Jk^  - k ^ , i f k 


Axial  slots 


dk  I 4>(m,k  ) F(m,k  )e“j^"^0 


, ab  cos(k  b/2)  |2 

■i  ( m . k ) = — — . z 

z 8R  (mij)  ) , ; 2 

(k  b/2)''  - (tt/2)^ 

^ I 

(2)  ^ 


5 


This  solution  is  suitable  for  numerical  calculation  if  (i)  < b for 

circumferential  slots,  and  < a for  axial  slots,  (ii)  kR  is  less  than 
20,  and  (iii)  the  medium  is  slightly  lossy  so  that  k has  a small 
(negative)  imaginary  part.  Based  on  this  solution,  extensive  numerical 
results  have  been  reported  by  Hughes  Aircraft  Company  at  Culver  City 
[7],  [8],  [9]. 


4.  EXACT  UI  MODAL  SOLUTION 


Under  the  one-mode  approximation,  another  exact  modal  solution  is 
given  in  [10].  This  solution  is  derived  from  the  Hughes  (SGI)  solution 
in  Section  3 by  a deformation  of  integration  contour  and  an  application 
of  the  Duncan  transform  [11].  THe  final  result  reads 


Circumferential  slots 


'>’,9  - G + jl 


^ cos  ni].. 


0 m--0 


cos  <f'.(n,k^)R(ni,k^)  dk,, 


cos  mil,, 

B = I 

m=0 


R(m,k  )G(m,k  ) sin  k z,,  dk 
z z z 0 z 


+ R(m,jn)ij)(m,jn)c  dn 

■'o 


2 k nik^^  j 

• r • ^ \irr-R  i 7--- 

L m t / 111  t 

= J^(X)  -f  Y^x) 

Iii  Til  m 


2 , m = 0 

1 , m 5^  0 


(4.5) 


(A. 6) 


) is  defined  in  (3.2)  and  k in  (3. A) 
z t 


Axial  slots 


r 


“ cos  m(J)_ 

0 Y U 

nkR  e 

m=0  m 


dk 


<1(01, k )e  —r, — - 

z 


'0 


N'"(k  R) 
m t 


+ j 


$(m,jn)e  '^^0 


dn 


where  l>(m,k  ) is  defined  in  (3.3) 
z 

This  solution  is  valid  only  if  > b for  circumferential  slots  and 

z > a for  axial  slots.  It  is  suitable  for  numerical  calculation  if 
0 

kR  is  less  than  20. 


5.  ASYMPTOTIC  SOLUTION 


(A./) 


The  two  modal  solutions  given  in  Sections  3 and  A are  based  on 
fields  in  the  Fourier  transform  domain.  An  alternative  calculation  of 
Y^2  involves  the  field  in  the  spatial  domain,  namely. 

Circumferential  slots 


‘12  ab 


dy^dz^ 


'^^2‘^^2  a^  yj^][cos  ~(Y2  ~ R'+>q)  ]gjj^(s,e)  (5.1) 


Axial  slots 


'"n  = i!  L ^^1^^^ 


dy^dz^  [cos  ^ z^][cos  - ZQ)]g^(s,0) 

2 

where  (y  , z ) = a typical  point  in  the  aperture  of  slot  n (n  = 1 or  2) 


n n 


= aperture  of  slot  n 


= Ay2  - Yj^)^  + (^2  ” ^1^^ 


0 = tan 


(z.,  - Zj)/(y2  - y^) 


(5.2) 


(5.3) 
(5. A) 

(5.5) 

(5.h) 
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Several  versions  of  the  Green's  functions  and  g have  been  approximately 

9 z 

determined  under  the  condition  that  kR  >>  1.  They  are  listed  as  follows: 
OSD  Asymptotic  solution  [12]  [13] 


g G(s)  [v(C)  sin  6 + cos  0] 


g^  -V  G(s)  ['/(C)  cos  6 + 


(5.7) 

(5.8) 


PINY  Asymptotic  solution  [9]  [14] 
g 'V  G(s)  r v(C)  [sin^e  + ^ 3 sin^e)] 


+ ^ sec‘^0[u(C)  - v^(C)  sin^0] 


(5.9) 


g G(s)  v(C)  [cos^0  i-  (2-3  cos^0)  ] 
z ks 


(5.10) 


UI  Asymptotic  solution  [ 15 ] 


g,  G(s) 


v(C)  [sin^0  + cos23]  + (r^)  u(C)[cos^0(l  - -^)  + (-^)  sin^0] 

ks  Ls  K.S  ICS 


ks  ks"  'ks 

2, ,-2/3.  . ...  . 2„  . 4„  . i . . ...  2 


+ 3 (v^  kR/cos  0)  [v*(5)  sin0  + (tan  0 + ^ ^ 


(5.11) 


g^  = G(s)  v(C)[cos^0  - ^ cos20]  + (^)u(C)  [sin^0  (1  - ^)  + (^)  cos^0 


ks" 


+ j(/2  kR/cos^0)  ^^^[v' (C)cos^0  + (1  + ^)u' (C)sin^0] 


ks"  'ks" 


where 


k Yq  -1 

= 2^ri^  ’ "o  = 

C = (k  cos^0/2R^)^^^s 


(5.12) 

(5.13) 

(5.14) 


The  Fock  functions,  u,  v,  etc.,  can  be  calculated  from  the  following 
two  representations: 

For  0 < C < 0 . 7 


8 


34“  -j3ir/4  1/2  . Zi  2 63/7  JH  - 2.A85  x 10  (5 


j3ir/4  1/2 


8.48849 


9.02265 


10.04017 


11.00852 


11.93602 


11.47506 


12.82878 


12.38479 


n 

1 

2.33811 

1.01879 

2 

4.08795 

3.24820 

3 

5.52056 

4.82010 

4 

6.78671 

6.16331 

5 

7.99413 

7.37218 



k 


r 

I 

f 

I 


It  has  been  verified  through  several  hundred  numerical  examples  that  the 
UI  asymptotic  solution  given  above  is  in  excellent  agreement  (within  a 
quarter  db  in  magnitude  and  a few  degrees  in  phase)  with  the  exact  model 
solution  for  all  slot  separations  ((^^,2^)  provided  that  kR  > 5. 

In  using  the  asymptotic  solutions  for  calculating  the  self  admittance 
care  must  be  exercised  in  avoiding  the  singularity  in  the  Green's 
function  which  occurs  at  s = 0.  A most  convenient  way  to  avoid  this 
apparent  difficulty  is  to  (i)  use  a large  number  of  points  for  the  two 
surface  integrals  in  (5.1)  and  (5.2),  and  (ii)  .shift  slightly  the 
integration  nets  for  this  two  surface  Integrals. 

6.  EXACT  PLANAR  SOLUTION 

In  the  limit  kR  -*■  “ the  Green's  function  of  the  UI  solution  in  (5.11) 
and  (5.12)  is  reduced  to 

g^  = G(s)[sin^e  + (2  - 3 sin^0)(l  - ^)  ] (6.1) 

g^  = G(s)[cos^e  + ^ (2  - 3 cos^0)(l  - ^)]  . (6.2) 

IJhen  (6.1)  and  (6.2)  are  used  in  (5.1)  and  (5.2),  we  obtain  the  exact 
solution  (under  the  "one-mode"  approximation  of  course)  for  two  slots 
on  an  infinitely  large,  conducting  plane. 

7 . APPROXIMATE  SOLUTION 

Based  on  the  UI  asymptotic  solution,  a simple  approximate  solution, 
is  reported  in  [10],  i.e.. 

Circumferential  slots 

Yi2  ' - [S(b  sin  0)  C(a  sin  0)]^  g . (7  1) 

7T  ^ 


10 


Axial  slots 


i 


[S(a  cos  0)  C(b  sin  Q)]^  g 

IT  ^ 


where 


s(x) " . CM  - 

1 - (kx/Tt)^ 


The  (simplified)  Green’s  funcLious  e.  and  e are  piven  by 

d)  ^ z ^ 


g = G(s)  [v(r)  [sln^  9 + t 


“ — COS  2G)  + u(C)  cos  G 

KS  ' KS 


kU  cos  0)  sin^  0 


v(^)  |cos^  ® ^ ^ °]- 


This  solution  gives  an  accurate  numerical  result  (within  several 
percent  in  magnitude  and  less  than  5°  in  phase)  provided  that 
kR  > 10  and  the  slot  separation  is  greater  than  two  wavelengths. 


8.  CONCLUDING  REMARKS 


(7.2) 

(7.1) 

(7.4) 

(7.5) 


Based  on  extensive  numerical  data,  we  conclude  that  (including 

as  a special  case)  can  be  best  calculated  by 

(i)  Hughes  modal  solution  if  kR  < 5 and  is  less  than  the  axial 
dimension  of  the  slot, 

(ii)  UI  modal  solution  if  kR  < 5 and  is  greater  than  the  axial 
dimension  of  the  slot,  and 

(lii)  Ul  asymptotic  solution  if  kR  > 5 for  all  slot  separations. 

If  several  percents  of  error  are  acceptable,  the  approximate  solution 
can  be  used  if  kR  > 10  and  the  slot  separation  is  greater  than  two 
wavelengths . 
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APPENDIX  A:  NUMERICAL  RESULTS 


DATA  SET  A OF  MUTUAL  ADMITTANCE 

(1)  The  mutual  admittance  between  two  circumferential  slots  on  an 
infinitely  long  cylinder  is  calculated  from  the 

* (Exact)  Hughes  modal  solution 

* (Exact)  U1  modal  solution 

* U1  asymptotic  solution 

* OSU  asymptotic  solution 

* PINY  asymptotic  solution. 

The  parameters  are 

* Frequency:  f = 9 GHz,  k = 4.7878  (inch)"^  A = 1.3123" 

* Cylinder:  R = 1.991"  unless  specified  otherwise 

* Slot  A:  Circumferential 

a = 0.9"  = 0.6858A 
b = 0.4"  = 0.3048A 

|y^J  = 1.70747  X io“^  mho  = -55.35  db 

Y = 1.8155  X lo"^  mho 
g 

* Center-to-center  distance  between  tvio  slots  is  (R(j)Q,Zy). 

(2)  Y^2  is  listed  in  (db  value,  phase  in  degree),  where 

db  value  = 20  logj^g  mho). 

(3)  Data  are  presented  in 


TABLE 

A-1: 

<))q  = 0 and  various  z^ 

A-2: 

Zg  = 2"  and  various 

1 

< 

Zq  = 0 and  various  ((>q 

> 

1 

(fy  = 0 and  various  Zy. 

Figure 

A-l: 

Mutual  admittance  Y^^2  betweeii 
as  a function  of  <t)Q. 

two  circumferential 

slots 

A-2 

Mutual  admittance  Y^^  between 
as  a function  of  Zy.“ 

two  circumferential 

slots 

A-3:  ^ cylinder  (UI  modal  solution)  and  that 

on  a plane  as  a function  of  z^. 

A-4 : Y^.j  on  a cylinder  as  a function  of  the  radius  R ot  the 

cylinder . 
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TABLE  A-1 


I 


Y,„  OF  SLOT  A FOR  (J)  = 0 

12  0 


■ 

Modal 

Asymptotic 

N 

Exact  Planar 

Z 

o 

Hughes 

UI 

UI 

OSU 

PINY 

R=oo 

1 

-62.62  db 

-62.62 

-62.54 

-64.22 

-61.7 

1 

-63.69  1 

0.3" 

-72° 

-72° 

-72° 

-43° 

-68° 

-67° 

-71.87 

-71.78 

-71.66 

-73.67 

-70.96 

-73.53 

2-. 

-117° 

-117° 

-116° 

-100° 

-118° 

106° 

1 

-82.3 

-81.84 

-81.83 

-85.46 

-80.80 

-85.4 

8" 

33° 

34° 

37° 

55° 

34° 

54° 

1 

-86.48 

-86.6 

-91.41 

-85.26 

-91.40 

16" 

-4° 

-1° 

20° 

-4° 

19°  j 

-91.95 

-92.46 

-99.34 

-90.83 

-99.33 

40" 

-115° 

-110° 

-83° 

-112° 

-83° 

b - - 


TABLE  A- 3 

OF  SLOT  A FOR  z = 0 
12  o 


TABL1-:  A- 4 


on  a cylinder  (UI  solution)  nnd  that 


VNVld  / ^'k  NIIAO ) 
®o 

I 


DATA  SET  B OF  MUTUAL  ADMITTANCE 

(1)  The  mutual  admittance  between  two  circumferential  slots  on  an 
infinitely  long  cylinder  is  calculated  from  the 

* (Exact)  UI  modal  solution 

* UI  asymptotic  solution. 

The  parameters  are 

* Frequency:  f = 9 GHz,  k = A. 787787  (inch)~\  a = 1.3123" 

* Cylinder:  U = 1.991"  , 3.777",  6" 

* Slot  B:  Circumferential 

a = 0.656168"  = 0.50X 
b = 0.01j123"  = O.OIX 

* Center-to-center  distance  between  two  slots  is  (R$„,z-). 

0 0 

(2)  is  listed  in  (db  value,  phase  in  degree),  where 

db  value  = 20  log^^  ^('^12^  mho). 

(3)  Data  are  presented  in 

TABLE  B-1:  = 0 and  various  Zq 

B-2:  Zq  = 2"  and  various 

B-3:  Zq  = 8"  and  various  4)^. 

B-A:  Comparison  of  Hughes  and  UI  solutions 
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TABLE  B-1 


UI  SOLUTIONS  OF  Y OF  SLOT  B FOR  (J 


Asymp  I Modal 


Modal 


Asymp 


Asymp 


-111.48 

36° 

3 

TABLE  B-2 


UI  SOLUTIONS  OF  Y, ^ OF  SLOT  B FOR  z 


Modal 


Modal 


Asymp 


Modal 

Asymp 

t 

' 

-106.89 

-106.94 

-177° 

-177° 

-115.80 

-115.93 

11° 

12° 

-124.77 

-124.95 

140° 

141° 

• 

-136.67 

-136.82 

106° 

105° 

-148.07 

-147.24 

51° 

44° 

-155.92 

-165.47 

-170° 

-102° 

Modal 


TABLE  B-3 

UI  SOLUTIONS  OF  Y,„  OF  SLOT  B FOR  z = 8" 
12  o 


Asyrap  Modal 


Asymp 


K 

Modal 


6.0" 

As  ymp 


-111.63  db  -111.7A  -113. A5  -113.47 


-114.44 


112.08 


-114.54 


-116.18 


-116.52 


-34 


34 


113.18  -115.94  I -116.26 


112.83 


-118.94 


-32 


130 


-129 


114.41  -115.12  -119.29  -119.82 


124.43 


116.70 


117.70  I -123.69  -124.22 


131.31 


27 

146.21 


-122.98 

. , ,0 


124.10  I -134.62  I -134.27 


TABLE  B-4 


COMPARISON  OF  HUGHES  AND  UI  SOLUTIONS 


Hughes 

Modal 


Hughes 

Modal 


Hughes 

Modal 


Modal 


Modal 


Asymp 


Modal  I Asvmn 


Asymp 


DATA  SET  C OF  MUTUAL  ADMITTANCE 


(I)  The  mutual  admittance  between  two  axial  slots  on  an  infinitely 
long  cylinder  is  calculated  from  the 

* (Exact)  UI  modil  solution 

* UI  asymptotic  solution 
The  parameters  are 

* Frequency:  f = 9 GHz,  k = A. 7877  (inch)~\  X = 1.3123" 

* Cylinder:  R = 1.991",  and  other  values 

* Slot  C:  Axial 

a = 0.4"  = 0.3048A 
b = 0.9"  = 0.6858X 

* Center-to-center  distance  between  two  slots  is  -^)  . 

^2)  listed  in  (db  value,  phase  in  degree),  where  db  value  = 2''' 


( |y^  , I in  mho) 


(3)  I'.i!  i are  presented  in 


TA3LE  C- 1 : = 0,  R = 1.991",  and  various  z^. 

C-2:  z.  = 1.5",  R = 1.991",  and  various 

0 0 . 

C- 1 • = 0,  z^  = 8",  and  various  R. 


Figure  C-1:  ^ cylinder  (UI  modal  solution)  and  that  or 

a plane  as  a function  of  z^. 


TABLE  C-2 
OF  SLOT  C FOR  z 

Modal 


DATA  SET  D OF  MUTUAL  ADMITT/\NCE 


1 


(1)  The  mutual  admittance  between  two  circumferential  slots  on  an 
Infinitely  long  cylinder  from  the 

* (Exact)  UI  modal  solution 

* Ul  asvmtotic  solution 
The  parameters  are 

* Cylinder:  R = 1' , 2\,  4X,  lOA,  “ (planar) 

* Slot  D:  Clrcumi  v.rential 

a = 0.3'* 
b = O.OIX 

* Center-to-center  distance  between  two  slots  is  (R4  ,z  ) 

o o 

(2)  is  listed  in  (db  value,  phase  in  degree),  where 

db  value  = 20  ^'^12^  mho) 

(3)  Data  are  presented  in  r- 


TABLE  D-1: 

= 

0,  R = 2>.  and  various 

n-2: 

S 

0 and  various  R and  z 

o 

D-3: 

\ 

= 

0 and  various  R and  z 

D-A: 

z 

0 and  various  R and  ({) 

o 

D-j. 

z 

=s 

IX  and  various  R and  i 

o 

D-6 : 

2 

o 

= 

5>  and  various  R and  i 

I 

i 
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TABLE  D-1 


UI  SOLUTIONS  OF  Y 


Modal 


Asvroptotic 


-110.84 

73° 

6A 

-112.19 

73° 

114.28 

-114.33 

2° 

73° 

j 

115.12 

-115.18  ! 

1° 

i 

-115.94 

IX 

j 

-98.60  db 

71° 

2A 

-103.87 

74° 

1 

3A 

-106.98 

75° 

TABLE  D-2 


r 


UI  ASYMPTOTIC  SOLUTIONS  OF  Y.  « OF  SLOT  D FOR  =•  0 

12 


r 

I 

'% 


TABLE  D-5 

UI  ASYMPTOTIC  SOLUTIONS  OF  OF  SLOT  D FOR  z = IX 


ii- 


TABLE  D-6 


DATA  SET  E OF  MUTUAL  ADMITTANCE 


^l)  The  mutual  admlttunre  between  two  circumferential  slots  on  an 
infinitely  long  cylinder  is  calculated  from  the 
* U1  asymptotic  solution 
The  parameters  are 

*Cy  Under:  K lA,  2A , 4A , lOA 

*Slot  E:  Circumferential 

a = 0.5A 
b = 0.2A 

*Center-to-center  distance  between  two  slots  Is  (Rdi  ,z  ) 

o o 


(.2)  listed  in  (db  value,  phase  in  degree),  where 


db  value  = 20  l^gjQ  ([y^^I  in  mho) 


(3)  Data  are  presented  in 

TABLE  E-1:  z = 

o 

E-2:  z = 


o,  various  and  R 
o 

0.5A,  various  and  R 
o 


E-3:  z = lA,  various 

o 


and  R 


E-4:  z 


2A , various  <P  and  R 
o 


E-5:  z = 4A,  various  and  R 

o o 


E-6:  z = R> , various 

o 


and  R 


E-7 : Comparison  of  UI  asymptotic  and  UI  modal  solutions 
E-8:  Comparison  of  UI  asymptotic  and  UI  modal  solutions 


E-1: 

Mutual  admittance  Y^,^ 
as  a function  of 

between 

t wo 

circumf  erential 

lots 

E-2: 

Mutual  admittance  Y,^ 
as  a function  of 

between 

t wo 

circumf erent ial 

slots 

E-3: 

Mutual  admittance  Y^^ 
as  a function  of 

between 

two 

c ireumferent ia i 

lots 

E-4 ; 

Mutual  admittance  Y^ 
as  a function  of 

bet  wei'ti 

t wo 

c ireumferent ia 1 

■let;. 

E-5: 

Mutual  admittance  Y^ 
as  a function  of 

between 

t\’o 

c i 1 v umf erent ial 

;.l,.t  s 

40 


TABLE  E-2 


TABLE  E-4 


UI  ASYMPTOTIC  SOLUTIO^B  OF  Y,„  FOR  z = 2X 

12  o 


o 

R = IX 

R = 4X 

R = lOX 

! 

-77.24  db 
70° 

-78.67 

76° 

-79.46 

80° 

-79.96 

82° 

10° 

1 

1 

1 

-77.52 

67° 

-79.44 

65° 

-81.77 

37° 

-89.55 

-148° 

20° 

I 

I 

1 

-78.37 

-81.60 

31° 

-87.38 

-79° 

0 

o 

-79.73 

42° 

-94.18 

112° 

-114.68 

-144° 

1 

45° 

-82.59 

9° 

-90.76 

-130° 

-128.88 

15° 

60° 

-86.17 

-35° 

-97.24 

94° 

-113.88 

177° 

-141.23 

153° 

44 


TABLE  E-7 


COMPARISON  OF  UI  ASYMPTOTIC  AND  UI  MODAL  SOLUTIONS 


Asym. 


72.54  db 


75.33 

55° 


86.37 

-77 


86.49 

76° 


103.77 
1° 


104.13 


TABLE  E-8 


COMPARISON  OF  UI  ASYMPTOTIC  AND  UI  MODAL  SOLUTIONS 


IX 

R = 2X 

Planar 

Asym. 

Modal 

Asym. 

(Exact) 

-61. 

-117° 

-68.69 

-114° 

-69.35 

-110° 

-72.28 

68° 

-73.64 

73° 

B 

-74.52 

79° 

-11 .Ik 
70° 

-78.98 

75° 

-78.67 

76° 

-80.29 

84° 

-82.10 

68° 

-84.3 

75° 

-84.01 

75° 

1 

-86.25 

87° 

-86.7 

66° 

-89.41 

72° 

-89.18 

73° 

-92.25 

89° 

48 


_l 

< 

Q 

UJ 

O 

(O 

< 

s 

X 

CL 

2 

• 

I— 

o 

>- 

X 

< 

O 

CO 

o 

X 

UJ 

< 

ir 

UJ 

M 

CL 

_j 

M 

CL 

M 

<I 

=) 

< 

=J 

s 

I 

1 

(T 

O 

1 

o 

Z 

1 

1 

Ur  EXACT  MODAL 


DATA  SET  F OF  MUTUAL  ADMITTANCE 

(1)  The  mutual  admittance  between  two  axial  slots  on  an 
infinitely  long  cylinder  is  calculated  from  the 
*UI  asymptotic  solution 
The  parameters  arc 

*Cylinder:  K = IX,  2X , 4X , 10^ 

*Slot  F:  Axial 

a = 0.2X 
b = 0 . 5 a 

*Center-to-center  distance  between  two  slots  is  (R*  , z ) 

o o 


is  listed  in 

(db  value. 

phase  in  degree) 

db  value  = 

20  lo:,^^.  (| 

\'i2  1 in  mho)  . 

are  presented  in 

TAiiLE  F-1: 

z 0,  varion^  ^ and  R 

• ) <) 

F-2: 

z - 0.  5X , 
r. 

various  * and  R 
o 

F-3: 

z "IX, 
o 

various  *.  and  P 
o 

F-4: 

z = 2X , 
o 

various  4i  and  R 
o 

F-5: 

z = 4X , 
o 

various  A and  R 
o 

F-6: 

7.  = 8X , 

c> 

various  j and  R 
0 

F-7: 

Comparison 

of  UI  asymptotic 

F-8: 

Comparison 

of  Ul  asymptotic 

F-0  ■ 

Compar i so- 

of  asvmntnt  ic  • 

Figure  F-1:  Mutual  admittance  Y^2  between  two  axial  slots 

as  a function  of 

F-2 : Mutual  admittance  Y^2  butween  two  axial  slots 

as  a function  of 

F-3:  Mutual  admiltance  Y^2  between  two  axial  r-lots 

as  a function  of  0^. 


5A 


TABLE  F-2 


TABLE  F-A 


TABLE  F-5 


UI  ASYMPTOTIC  SOLUTIONS  OF  Y^^ 


FOR  z 

o 


4X 


r-H 

11 

0^ 

R = 2X 

R = 4X 

R = lOx 

0° 

-111.54 

-178° 

-111.53 

-178° 

-111.52 

-178° 

10° 

-111.78 

177° 

-111.97 

168° 

-111.81 

132° 

-105.41 

-27° 

— 

20° 

-112.38 

164° 

-112.36 

126° 

-108.23 

21° 

-100.03 

-35° 

30» 

_ . 1 

-113.06 

143° 

-111.16 

67° 

-104.80 

-104° 

-100.65 

-154° 

B 

-113.38 

97° 

-108.48 

-24° 

-103.49 

28° 

-105.31 

100° 

60° 

-112.64 

47° 

-107.29 

-123° 

-104.94 

114° 

-111.46 

-67° 

59 


(ti 

R = 

1 

R = 2X  ! 

o 

Mod  a 1 

Asym. 

Moda  1 

Asyrn.  j 

0° 

-87.06  db 
-171° 

-86.65 

-173° 

-86.83 

-172° 

-OD . d3 
-172° 

10° 

-87.69 

176° 

-87.3. 

172° 

-88.23 

IjV 

-87.92 

134° 

20° 

-88.91 

139° 

-88.37 

128^ 

-87.64 
3 5° 

I 

-86.^1  1 

26  i 

1 

30° 

-89.40 

85° 

-88.07 

70° 

-87.01 

-72° 

-85.77  ! 

-81° 

45° 

-89.19 

-87.32 

-15^ 

-88.67 

119° 

1 

-87.30 

109° 

60° 

-89.84 

-83° 

-87.72 

-99° 

-91.86 

-61° 

-90.36 

-72° 

TABLE  F-8 


COMPARISON  OF  UI  ASYMPTOTIC  AND  UI  MODAL  SOLUTIONS 


R = lA 


R = 2X 


Planar  ' 


Modal 


Asym. 


Modal 


Asym.  (Exact) 


TABLE  F-9 


COMPARISON  OF  ASYMPTOTIC  SOLUTIONS 


PINY 


UI  Asynw  PINY 


Fit;urf  K-2:  Mut.i.il  iidmittance  Y , hi- tween  two 


Ittance  Y,„  between  two  axial  slots 


APPENDIX  B:  COMPUTER  PROGRAM  LISTING 


I 


This  appendix  contains  the  program  listing  of  all  solutions, 
except  the  exact  Hughes  modal  solution,  discussed  in  the  text. 


onrvo  OoorionrioonnnoODnorinonnnoonrjoonooriortonrjoriOonooononnn 


! 


ASYMI'lO'ilC  SOUiTIONS  OK  Y, 


Turs  i’ROCfi)i,rt  1$  ustp  TO  ?OLVE  iro/(  MUTUAL  ADMITTANCE  OT  DOTS  £lTUlR  OK 
A CiLlHOSA  OR  ON  A t^ANF 

TUERe  ARE  TWO  T^l’Es'OK  SlorS:  1 T 0^  8 CU  MF  E R EHT  I M ; 2)  AX  i At 

TMIi  V'PCGRAM  (NV0LVC5  A LOT  Of  I Mf  RC  R A T JC  K S WHICH  A^k  D'SICALLY  SOLV(-0 
t)T  5uA.MATtON  IIETHOO 

« « 

' r.V(»Ur  pAflAHeTERS  c this  program  * 

• • 

IPL,4!i-CO,17R0L  -THE  PROGRAM  III  DEALING  WITH  2 Um'FREI.'T  CASES; 

1)  IN  PLANAR  CASE,  SET  IPLAN=] 

2)  XA’  CYLINORICAL  CASE,  SET  lPLAN-2 

eOLLOvlING  PARAMETERS  ARE  CCMMON  TO  BOTH  IPLAN=1  AND  JPLAH-2 

CUK  -ASSIGN  A LOGIC  VALUE  'TRUE'  IF  HE  APE  IKTERFSTED  JN  THE 
CIPCUHFERENTI  AL  CASE;  OTHIPWISE  ASSIGN  TAISE'  TO  IT 
AXIAL-ASSIGN  A LOGIC  VMUE  ’TRUE’  IF  WE  ARE  ISITERESIED  IN  THE  AXIAL  CASE; 
OTtlfiRWlSe, ASSIGN  ’FALSE’  TO  IT 
ft  -THE  TONCeP.  l.EIIGrU  OP  A SLOT  (MEASURED  IN  'WAVELENGTH) 
a -THE  SHORTER  LENGTH  OF  A SLOT  (MEASl/REp  IN  KA  V RLE  HGTH) 

Iph  -number  of  SUH divisions  OVER  THE  LCNGtR  LENGTH 
IPa  -NUMBER  OF  SUBDIVISIONS  OVER  THE  SHORFfcR  LENGTH 
Ul  -01  = 1 IF  KB  use  Ut  ASYNP’i  Or  IC.  IF  NOT,  UJ=0 
D5J  -0SU=2  IF  WE  USE  OSU  ASYMPTOTIC.  IF  WOT.  OSII  = 0 
PINT  -PINY=)  IF  we  USE  PINY  ftSYMP/fiTIC.  IF  N6  T , prNY  = 0 

7.n  -Self  adhittahcc  of  Slot  (measured  ih  mho) 

Z -ARRAY  OP  20  elements  AT  MOST.  EACH  ELEMENT  STAUDS  FOR  THE 

SEPARATION  BEIHEEN  2 POINTS  ALONG  Z AXIS  (MEASURED  IN  WAVELENGTH) 
-NUnaCR  OF  ELEMENTS  IN  . DON'T  BE  GRt'ATEK  THAH  '/O 


FOLLO-WINC  parameters  are  only  good  for  CYCIIUDRICAL  CASE. 

jusr  foRuer  them  if  we  are  dealthc,  with  The  planar  case 

RADIUS  - array  of  different  RADII  OF  A CYLINDER  (MEASURED  IN  WAVELENGTH) 

NOR  - NUMBER  OF  ELEMENTS  TN  RADIUS.  DON’T  3E  GREATER  THAN  6 

PHI  - ARRAY  OF  6 GLFHENTS  AT  MOST.  EACH  ELEMCNT  REPRESENTS  THE  ANGULAR 

SEPAHATtON  BETWEEN  TWO  SLOTS  (MEASURED  IN  DEGREE) 

NO?Ht  - NUMBER  OF  FXEHENTS  IN  PHI.  DON'T  BE  GREATER  THAN  6 

rOUOVING  PARAiMETtlRS  ABE  FOR  THE  PLANAR  CASE 

YP  -ARRAY  OF  20  eLfiMENTS  AT  MOST.  EACH  OF  THEM  IS  THE  DISTANCE 
BErweSN  2 SLOTS  ALONG  Y-ftXTS  (MEASURED  IN  W A VEL  ENGT II) 

;/DY  -NUH0ER  OP  ELEMENTS  IN  YP.  DON ' f BE  GREATER  THAN  20 

I/1PLICIT  COMPLEX*  16  (C,H,Z) 

IMPLICTT  REAL-e  (A- H , D-G, P- Y) 

real  PHI  (6)  ,RADu1S(6)  ,Z(iO)  .FREQ (6)  , YP  (2D) 

INTEGER  PtNY, OSU, Uf, test  ^ 

LOGICAL  tOrt.AXXAL 
(\EflL*9  KA.MAG.Z* 

Real  inch 

mavAvmiv^^  0 , n 

CoMrtCN/DRTA 1/TN, TNpT, RRO, C 1 , C2 , F2, lOP 
COMMON  PI.TZI^TZi ,tV  J ,TY2. R#THEThA 
COMMON/>ATA/AQ.  BO, zo, yo 

BEFOfte  EACH  PUN, CHECK  THE  FOLLOHING  PARAMETERS  A^D  HAKE  APPROPRIATE  CORft.E.CTfON 
SET  ft  vAuie  TO  rue  kormaltzatiom  factor 

if1=1. 
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Ow'.  i(  4 iOw(  n w I iU.  ii  ..  j.  ji  i(  i(  J<  j(  H it  i(  i.  i ^ t 


onO 


C the  ^p?ftOPRI^Te  LOGIC  VAI-'lfS  TO  A^5AL  A 0 CUK> 

CU«^.  fALSe.^ 

c ser  A*THE 'length  or  longer  side  of  slot  and  b-tme  length  of  shorter  one 
A -0.5 
8-0.2 

c choose  a proper  intercration  grid 

I PA=1 « 

c choose^which  as  VHP  rone  he  trod  to  3e  used 
1)1-0 
050=2  ^ 

P I N V = 3 

C assign  the  NU/10EP  OF  DIFFERENT  VALUES  OF  StPARTION  8EIWEEN  2 SLOTS 
C ALONG  2-piRCCTCON. 

C then  construct  a corresponding  array  Z of  J<0Z  ELtHLNTS 
N0Z=1 
2 ( 1\  = 0. 6 

C SET  I^LAfi=1  rOR  PLAKAR  CASE  AND  IPLAN  = 2 FOR  CYLINDRICAL  CASE. 

c IP  iPlah^i.  skip  this  section  to  the  next  OHE 

C IP  lPLAN-2!  SET  THE  NUA1BER  Of  DIFFERENT  RADII  OF  CYLINDER  TO  NPf? 

c AND  then  Construct  a corresporoivg  array  radius  cf  ndr  elenents 

c SET  NUMBER  OF  DIFFERENT  VALUES  OF  ANGULAR  SEPAFATIOK  atTHEEN  2 SLOTS 
C THEN  CCNSIRUCT  A ARRAY  OF  PHI  OF  NDPHl  El FMtNTS  (MAX.  BO.  OP  ELFMEKTS  IS  6) 

raSiusM) 

RADIUS  2)  =2. 

RADIUS  3)  =4  . 

RADIUS  4) =10. 

NOPP.I=^ 
p H r ( 1 =0 . 

PHI  2 =10. 

PHI  3 =20. 

?HI m'  = 30. 

Phi  5 =45. 

PHI  6 =&0. 

C IP  tPLAK=i  YOU  CAN  PUT  THE  WHOLE  DECK  INTO  THE  READER  NOW 

C IP  IPIAN=1.  SET  THE  NUMBER  OP  DIFFERENT  VALUES  OF  DISTANCE  EETHEEH  2 SLOTS 
C AlOHG  Y-OIKECIION  „ . , . „ 

C tHEj  CCNSTRUCr  A CORRESPONDING  ARRAY  YP  OF  NDY  LLEMENTS 

NDY=1 

Y P (1)  =0.  000  1 

APIER  HAXlNC  THE  CORRECTIONS,  YOU  CAW  PUT  THE  DECK  INTO  THE  READER 

CCUM=1 
CAXIAL  =2 

tP(.NOT.CUH)  ICUH=2 
IF(. NOT.  AXIAL)  IAXIAL=1 
CL=1 

TEST-PJNY-OSU-Ut 
CP  (TEST.  EQ.  3)  IL=3 
CF(TE5T.EQ. 1.0R.TEST.EQ.-2)  IL=2 
IU  = «AXO(PINY,OSI1,  UI) 

P I=4‘ OATAN (1 . DO) 


GO  TO  3 
vJRITE  (6,666) 
xIPITC  (6.  S55) 

^K  = Z.  *PI 

DO  3)  lU  = ICUf1,  ^^)(IAL 

S8-I2,P 

8o=e»x5 

GO  TO  6 
5AVE=  A 

VIavc 
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\ssi.  ;•)  "II'-'  \ t’.M  i.o(;:r  v't.'i;.'  to  av:,’'!  At.n  c".'-- 

:i!  = . bM.s 
\ T T A L . ■’ '!  I ; . 

3n:  1.  jp  mj'.cp-  op  piot  avu  r”''  .s'lrPTK'’  o:i 

. ■' 

!i  - ^ . P 

' ijj  : a"p'  j;*-;«  ;:ir  ;p 
A ■ 1 !4 


p.nj-i-:  A.i::i  -v 'y  i ' ’’p::"  '"■-•kji)  "o  '’i:  uppp 
;:  = ■■) 

J j!I  -p 
P IN 'A  - 1 


i ••  iJii;’’  O'-’  v.’-in;-;  o»'  :svi-r, r p.lo'"s 

r-ij.  TT'::'"iTi'  ' co='-'i  i-opornA  no  a')pm  z or  i d?  m.in-'pnts 
r)"-=i 

.$■-’.1'  -'1  Pl-\n\i:  CAN’''  ano  !>■’!,;, n-2  ro"-'  ryiiri'P’c.’i.  cast. 

: > !,  A ; - ■’ 

•- p >i  • - T r'*  " ) tut;-  ON'-' 

' ' I’A'i-"'  N " It 'N'MP;"'*r  p p*  ryi  to  n^'P 

\ .'.)  ■•>  o'lpr .-p ” pi  ni;:o  ?'’:-ay  rv  pliti'-'nt’p 

; ik  ^ t."'’  ANon’^Ar  ppp'-pa"'"-^'-: 


A r:..  'Ti  'Jc 
))>=-! 

'j:  I : ( I) 

< "I : I N I 


0"-  "I’X  OF 


("V'.  NO.  OF  !;i.kmf;;tf  : 


<M)1  I ? {M  =■*. 

; (i)  =■'  \ 


p i:  m 
p i:  »'i  = ’o. 


’ti:  =-<-j. 


p ii  . 


Ir’  .P;. Y3I  PI"'  TH:-'  IK  T.  F ri"CK  FNTO  THr  rFAn'-"  liO'i 


I-T  A.iK  'MMP.FF  fV'’  DTFFiFi.NT  VM'IF: 


A Y- j" " :i 

r.iN,.  c;- 1 1 -'’r-'ic:  a L'OF?":'rr!;.r)‘"Nr;  apfay  yf  or  ’’dy  m.:- i';’:n  f 


r,IF'"’'-NC:-'  FPTV'l'HN  2 


iOY=  1 
f P (1)  0 jAI 


,\''riP  anz  co pp.f.:''’tovf,  yon  can  pm'"  ti'f  r’'icf-'  into  the  fkade;! 


'■AU-1^1 

-2 

: ? ( . 'I  j~.  c n)  :cn  -i=2 

I'-’  .'I  ax:;.-.,)  ia.xiai.=i 


A'-^FT-PA'IY-  )J  I-HT 
: ('A  -.' -^T.  A,).  ’)  IL-? 

IF  ( A.-:  I-.  F^.  1 , "'FFT.  r 
PSP,  II’-) 


0.-'’)  IJ.=  2 


P:=ii  i>  OA  ’A  I ( 1 . NO) 

; p tj  (1,  ?)  ,:i>la“i 

INIA-J  (N,  /P'’) 

;o  -"0  ’ 
i T A - I ( A , ''i ')  A ) 

-i  (AT ' -;  ( '■  , 5 A 0 5 
.•:r<=p.  *;■: 


))  ))  ri"i, '•« 'PiA  i 


I 


LO^LiJZ)  ’■KK 
C P ( I P \AH  . £Q  . 1 ; 


MO  PHT^ 


KO-K''X\< 

OO  ■ ir  k'K 
[rrr.p  -iPO 
L Pi^  ^ r PA 

tPA  =TEE11P 
W lOTh  ^-^AO/IPA 
v*i6rH2-30/rP3 

CI^COKyp  (OC/lPl.X  {O.0P,-l>T/3.)  ) 

C2=COPXP(DCPPLX  (0. 00,  Pl/<» . ) ) 

P2^D5QRr(?I) 

Y 1 ^-AO/i. -).'10TH  1 /2. 

21=-00/2.  -\WIO TH2/2. 

GO  rO  (7,3)  , IPLAN 
00  UO  lU 

DO  50  I«A0^1,flDR 
PPO  XK«/^A0I  JS  (I  RAD) 
h’RI  IC  (6  , U<4U)  XK 
IPfrj.PJ.2)  GO  TO  9 
WAfie  (6,33  3) 

GO  TO  11 
write  (6, 2'?2) 

WRCTE  (6,111)  A , AO  ,8,00 
IP  (IPLAN.  t’Q.  1)  GO  TO  12 
WRire(6,99)  ftADIUS(IRAD)  ,RHO 
GO  TO  (9  1 ,92,9  3)  , TOP 
WRITE  (6,1233) 

GO  TO  12 
WRITE  (6.12  1U) 

GO  TO  12 
w ft  r re  (6  , 1 2 3 5) 

WRITE  (6,  6 6)  I PA,  I PD 
DO  60  J2-1,HDZ 

CP  (ll^LAi.  EQ.  1)  MOPHI^NDY 
OO  70  C(=1,NDPHI 
CP(r?LAM .RQ. 1)  GO  TO  13 
YO-RHO 'P I ‘PHI  (lY) /ISO. 

IP  { PH  £([().  EO.  0. 0)  YO^O.001 

Y^YO/KK 

GO  TO  ia 

Y^-YP(IV) 

Y0=  Y»  XK 

WRITE  (6,77)  PHI  (T  Y)  ,Y,Z  (.JZ)  , ZO 
Y2  = Y0  -a6/2.  -Vi  IDTH  1/2. 

Z2  = ZO-80/2.  -wcr)rH2/2. 

ZSUfI=0. 

00  60  K=1,IP\ 

TY1  =Y 1 +H16tH1 *K 
00  90  L^1,IP3 
rii  = z 1 fvi  rC)TH2'^L 
DO  100  H=1, IPA 
TY2-Y  2iwrorHl 
DO  1 1£>  N = 1 , IPS 
T22^Z2tWIDrn2  >‘N 

ft^DSQftT]  (TY2-TY1)  >'>‘2+(TZ2-TZl)''*2) 

THETHA^DATAKi  ( ( TZ 2- TZ 1 ) t ( T ^ 2 - TY  1 ) ) 

1F(  IPIAN. £Q.  1)  GO  TO  6Q^ 

CALL  CYLIND  (IJ»  ZSUM) 

GO  TO  110 

CALL  PLANAR (IJ, ZSUM) 

C0HTIf4UE 
CONTINUE 
CON  rlNUE 
CONTINUE 

iY2-Z  SUM  *■  6N  [OTH  I *wr  0T|12)  +*  2'*  (-2. ) / ( A0*30) 
MAG--CDABS  (ZY2) 

PHASE^OATAM  2(01  HAG  (ZY21  ,OREAl  (ZY2)  ) <‘190./PI 
Z EX  PON^COC  AP/OCITPUf  (0.  oO  , 05Q.RT  ( YO’'''*  2 + ZO**  2)  ) ) 
ZPROO=ZY2'ZE<POM 

PHN^0ATMJ2  (OIMAGZ  (ZPftOD)  .OREALZ  (ZPROD)  1 *180. /PI 


‘*2*(TZ2-TZ1)**2) 


I) 


rilon 


GO 

50 

•tO 

JO 

66 


77 

73 

1 »1 


coil 

t r ( £p  LAN  . 1 ) 

COWrt  HUE 
COLTCNtlE 
CONTINUE 
fORHATC  ' 

1 ilWlSE)  V 
i 'O' 

L’ORKAT  ( 

FOr.'JATdx 

PQft/IAT  (' 

1 JXj 


i,;o  TO  30 


)X,'*  NOnrtAI.  I /ATION  : AOS 


})(,'•  integratton  cpid 
ox,  '♦<■«*♦  DATA  CUrPlJT  * 
?U1=  


(UNl ESS 
X',r3/ 


SPECIFIED  OTIIE 


'Q','?UI=  ',E0.7,  • <DEG>:  ' ,5X,  ' y=  *,D15.7,'  <«AVEI.E 
'.ElO.u,'  «'WA  VFLrNCfn>;  • ,Sy,‘ KZ=  '.E10.N) 
}X,'rI2=  <HnO>',F7.?  » <DEG>; », SX, '0B= 

';  norm  PHASg=  >,n.2,  ' <bEC>') 


ENGTH'>  ; 
E12.5, 


t'Or.MATC 


«KR  = 


2?2 
3 3J 

U(l>4 

S5S 

ObO 

777 

oed 
1 > 33 
'2  3i 
12  35 


POPfiA  T ( 
PORMAT I 
FOR  NAT 
FOR  HA  T 
FORMAT 
PORfAT 
FORAU  Ti 
FOR.H  AT 
format  I 
FOPHAT 
STOP 
Cno 


CYLINDER:  •,5X,'P=  '.EIC'.H,'  < W A VE  LE  NG  T l!>  » 

£10. u/'  ',3X,'»  METHOD  OF  SOLOriO:!:') 

510T  DirE.NSrov:  •/ 

' i 6X,  • A= ',  E 10.  u , • <VA  V ELENGTH>' ,5X,  • KA=  •,P.10.a/ 

' ,6X,  ' 0=' , e 10.  U , ' <HAVfU:NC.TH>  ' , 5X, 'KB=  ',E10.ii) 
AXIAL 'J 

*,3X,'*  CIRCUMFERFNTIAL') 

■ '■  ■■  K=  <1/W.T  VEIEHGTH>' ) 


' - • ?X  ' ♦ 

I Q I ' lO^l  l*^*'»**N 

'0','lOx!*''  MUTUAL  admittance  OF  A SLOT 
'O'.IOX,'”  mutual  ADKITTAHCE  OF  A SLO 
1OF0. 5/)0Pe.  5) 

' t'  ,Z6K.  ' UT  ASYMPTOTIC! 

'♦',26X,'  OSU  asymptotic') 

' »' ,26X, 'PlNV  ASYMPTOTIC) 


*V0  000  ft  ftftftft  *00000*0  0000**00  1% 

ON  A CYLINDER  *’ ) 
' CH  A PLANE  *') 


THIS  Su3R0’JriNE  IS  ysEo  ro  calculate  the  functions  CVF,CUF,CV1F,CVPF,CUPP 

S U3R0UT  IMG  FOCK  (X) 

IMPLICIT  COrlPLE»'16  (C,Z) 
implicit  ReAU-3  (A-R,0-H,P-Y) 
rem.  tk  (10)  ,tnp:  ( 10) 

C0MPLEX«16  dchplx 

Connori/CF/c  vf,CUF  ,cv  1 f‘,cvPF‘,  cup? 

COMHON  pi, TLI , TZ2 ,TY1 .TY2, P ,THETHA 
COMMCN/OATAl/rN,'t'NPI,KUO,C1  ,C2,  F2 

k0nm., 

C VF-J . 

CUF^O . 

CVIF^O. 

svpr-0. 

CuPfrO. 

PO  io  N^1,10 
'Cl 

ZXNPl=TH3l (N)  *C1 
C 3 = CD  £Xp  (DCnPlX 
CU-COEY.P  joCKPLX 
C VF-CVFKTJ/ZrNP 
COf  ^CUF»C>4 
CV  1 F-CV1  fi-CT 
cvpp=  ( 1. DD-DCKPLX 
CUPF-  (1  . DO-OCMPLX 

COHtt  HUE 


<0, 

r 


00, -X)  *1TNPI) 

00, -K  *arN) 


2<J 


C V p F 

(iciuii 


I HUE 


(0. 

io. 


D0,2*X)  * 

DO,  2 *X/3 


ZTNPI)  TC3/ZTNPI  «-CVPP 
.)  -ZTH)  *Cu  + CUPF 


CVF-Fl’F 1 'CYF/C2 
CUF-1. ' FZ'f  J‘C2*CUF 
CV1  F^  Zj.  «(i'F3  »C2*CV 


CvPF/(2, ‘FI 
F2 'FI “CZ*CUPF 


ESTRY  fOCVl  (X) 

F 1 = t>SCjftT  (X) 

F 3 X * * 3 

2l--F2*CZ»0S(^RT(fT) 


71 


non  non 


:2  DC.'IPLX  (0.  DO,  1 . 00/0  0. ) '3 

13  --  * h.  /2. ) / ) 

F S ■ F i • ' 2 

L'V  K - 1 . D0-/:1  /u  . ♦ ]<■  Z2  *7.  * '/.l/d,  "i.  1 HI  D-S^FI/ 

ClIF^I  .D0-F.1/’.  ♦2>.‘7-2»5.  *■  ’ ■’  ' 

CVIr^  ’.OOtZl/2.  -1^).  '■'ll-  '. 

CVPFO.  *F2«2V(3.  *C2*"31  *■ 

C 'J  P F = 3 . • K 2 * K ' / ( H . 'C  r » '<  3 ) ♦ 

( 2210-1  >F«/X 

RETURN 
END 

THIS  jIlOROUTCNE  IS  US£0  TO  GET  THE  'PLZNAR*  SOLUTION 


l,  ).  ! 'J  ! lI-  t 

*Z3  «'i.  '■;55D-2*KH 
71.  ’7  2/X*63‘zVl16 
OC/i  PlX  (0. 1'0,5.  ■>  /'■ 


6.  •'XI 

ocj 


-2. H 85D-2*  FH/x 

*23/  (2*x) 


n 


SUBROUTINE  PLAN  AR  (I.I.ZSUH) 

ITlPLl  Cl  r COaH>LE X-  ’6  (C,  H,Z) 

IMPLICIT  REAL'=  3 (A-0/ 0-G,  P- V) 

3 E ft  L * 3 Z 3 

COI'.HON  Ft  ,TZ!  ,T22,TY1,TY2,  R.THETHA 
C0T.M-).l/0ATR/A,n,T5,  YO 
GO  ro  flo,  30)  ,I  J 
XKl  = ( (TZ  2- TZl)  /R)  •*  ’ 

Xh2=2.-3. *<1l 
XPL  = X 1 V3 
XIH-XM1*X'13/R*»2 
Hrt=Cf)tXP  {OCNPLX  /0,DO,-P))*DCHPLX(XRI..-XIf' 
/ACTORS PCOS (?I * :Y 1/ft) *DCOE (PI* (TY2-Y0)/A) 
1 3 U K ^ ZS  UTl » F A C T 0 H ' H A 
.1ETURN 


-XIF) /( 2H0. * R*PI*  *2) 


?0 


ZAl=0CHPLX(1./3**2,’./R1*(?.-:- 
ZA2  = C HEX P (DCrtPLX (Q. pO, :P)  ) * ( DCO S (TMFTHA) 


HA=  (0. ,-1.)  *2A2/  2U0. 
?ACTOR-=OC03  (P-‘TE  V0) 


“DCOS  (THETEA)  ‘*2) 

A ..  *.2*2A1)/P 


Z3UR=  ZSU.RtFAI 
RETURN 
END 


rOR»HA 


PI ’'  *21 
ncos (PI* (TZ2-Z0) /6) 


THI3  SUBROUTINE  IS  USED  TO  GET  THE  'CYLINOEICAl'  SOLUTION 


10 


40 
2 1 


SUBROUTIIIE  CYLIND  (I J,Z5fIM) 

IMPLICIT  C3MPLr.X‘16  (C,H,7) 

IMPlICir  REAL'S  (A-8, D-G, P-Y) 

REAL'S  ZO.XA 

REAL  TN  ( 10)  .TtlPJ  (10) 

COM  HON  P:,T21 ,Tzi ,T Y 1 , T Y 2 . P , THFTK A 
COHHON/CF/CVF,CUr'  ,C  V ) F,  C V P F , CU  P '' 

COM  HON/D  AT  A 1 /TN  , T HPI  , RliO  , C , C?  , F2,  lOP 

COHMON/|)ATA//i  ,P,iO,  Y<5 

ZGR=  (0.  1.  ) *CDfX  P(DCHPLX  (C.  PO.-P)  1/ (2U0.  *P!-PI**2) 

AHGLE^DATAN2(0ABS  (TZ2-T/.1)  .DABS  (Ty5-TY1))  *180. /PI 
I7(ANCLC. tr.33. 99)  GO  TO  lO 
I T=1 

THE THA=PI '39. 99/1 RO. 

3H=  ( 1 . , 1 .)  /pFCDEX  P(DCHPI.X  (0  . DC  , - PI/9  . DC  ) ) *DSOKT  (PI  *R/2  .)  /RllO 

RHOG--RHO/DCOS  (THETHA)  **2 

KA  = P*  DftnSJ  ( 1.  /(2.  ♦PHCG*»2)  ) **  (1  ./3.  ) ) 

IF  (KA  . I.T.().  T)  GO  TO  20 
CALL  fOCA(KA) 

GO  TO  30 


CALL  F0CK1 (K\) 

CP  (IT.  EO.  I)  GO  TO  90 

(Q.  , 1 . ) • (1  ./DCOS  (THETHA)  **2J  /P*  (CUF-C  VI  F"  DSlK  (THETHA) 
CO  JO  (21,  22, 2 3)  , lOP 
ZB81-  DCrtpl.;^  ( I.DD,  -1.  /R)  *CVT 
2llf2  = CUr 
2hl 


2) 


ZhR9*  JO.  ,’.i  ' **  ' cuPf*i)TM.'(THFTHA) 

HB-ZGR* (Z m -*H92«ZHy3*ZHFU) 

2TM1*  (0. , 1. ) *CU  PF/(DSQRT (2. DO) * RHOG)  •» j 2. /3. ) 

HT  = (0.  . 1 .1  * IGR/R*  (CVF'DCMPLX  ( 1 . 00»  - 2 ./!  ) ’I'UF4ZTA\1 ) 


r 


- >-:a  *!>cos  ( rH£T;.’A ) • • 2*)tr*nsi  r;  (1  nirHA)  * • 2 

ilVH  ill'  £>5^H  (TH  tT'lA)  *-2*HT‘  OCOS  (Tl! f T Hf.)  • ‘ 2 
G'j  ro  500 
iG  A"cvi' 

rtT  (0 . . ’ • ) "-cj r/p 

HZ  \b*0cO5  fTHETHA)  * * 2 ♦ HT' 05 1 M (T  HE  Ti;  A)  < » 2 
tiPH  t MIB  ' 0$c;f  ( TliETHA)  • *2*  iiT*  5CC-;  (li)FTHA)  • »2 
;o  TO  500 

Z2H2  = - 3-  * n^-lETHA)  2)  /R 

;l  2 = 7.G  rt'C/E  ' (OCO  3 ( TMETP.A)  'T-*  (0.  .I.  ) * (?.  -3. 
Hl’H  r ■■  ZCfi*  (C  VE  » ( DSi:i  (TiltTMA)  * *7*  i JH?)  «ZW) 
50)  ZGKECt;=  IPMt 

IffIJ.EQ.J)  S3Pl-eK=HZ 

FAC  ro PCD  5 (pi'-r r i/a)  ’ ocof.  (pi-  ( rv?- yO)  /h) 
tF  (IJ  . EO.  2)  FAC  rOP='j;05  (FT'  VZ  1/B)  * DCOS  (Pi< 
zs'J.T -zoari*  'AC  i of'^  zgrf.cV 
RClijS  •< 

St.'D 


22 


2 ) 


DCOS  (7H  tniA)  ♦ * 2)  /R) 


{TZ2-7.C)  /Irj) 


I'l  I XACr  SOI-UJ  IUN  OF  Y 


12 


J i’ROO«f>K  fOK  CO  MPLirATlON  OF  THE  fiUTUAi  A DN I 'i  T A NC  E BETWEEN  IWO 

C li».,iriCAL  Clf!CU.‘ir£Kf;iri  AL  SLOTS  ON  A C Y L I K I;  £ f(  ( UI  MODAL  SOLUTIO;.') 

ICAL  KO,  K Z,  ^r  , I 2 , K/.K  I PC 
CvOm  PL  EX  » 1 b 2 1 , Y 1 2,  t SI  f.X  P , Y !)1  2 

e'  I (400)  , ?K  (<.C0)  , FN  (4CC)  , DIMAG, DSEAL,  DATAN2 
C IM'JC  PASAiUrESS  ; 

C KO-W.WE  i.  J'MSEd  IN  FB££  SPACE  IN  lEilMS  OF  1/INCH 

C ■|'liO  = .lADIJS  JP  CYLINCER  <INCH> 

J .t''3=  SL3I  DIMENSION  , A>B  <INCH> 

C Pit  l0-=  ANGULAR  SEPASATION  OF  THE  SLOTS  (CENTER  TO  CENTER)  <KADIAK> 

0 ZO=  SEPAtiATiON  CF  THE  SIOIS  IN  2-DlrECTICN  <INCH> 

C ni=  NCRMALISaT  ION  FACTOR 

0 tUAX=  MAXirtUi  NC.  OF  TERMS  WHICH  HAS  BEEN  USED  IN  CALCULATION  OF 

J INFiNIIE  SERIES 

0 NCYCLE^NO.  of  subsections  QETWETN  any  two  SUCCESSIVE  ZEROS  OF  INTEGRAND 

C IM  rRfiPESoTDAL  RULE  FOR  NUMERICAL  INTEGRATION 

Pl= J.  141  59203 
'Y0=  1./  (120.  *PI) 

SO-2. *Ft 

FHcQO.310»KJ/(2.  »PI*2.54) 

.1  =0.3 
!1  = 0. 2 
AKa  -KO»A 

3 K D ^ K 0 • 3 
MMaX^  330 
hCYCLE=40 
Y 11  = 1 . 

RM0=2. 

fiK=KO*RHC 

v*  Ri  1(1  (6  , JO) 

3i  FOF.HaT  (401V,  ’MUrUAI  ADMITTANCE  OF  SLOTS  CK  A CYLINDER'/) 
w R 1 re  (6, 5 0) 

30  FORM!  T (52VC,  'ClftCUMFf  RENT  lAL'/) 
v.RirJ(6,5lj  FREC.SO 

5 I FORMA  f (2  5a,  'FREQatNCY ', 10X,'F=',E12,5,'<HZ>',10X,'K=',E12.5, 

^ • <1/INC»U'/) 

lYd  I re  (6,52) 

52  FORMAr(l0x, 'SLOT  DIMENSIONS') 

WRriC(6,53)  A,AKA,B,3K3 

J3  FORMAT  (1  OX,  'A  = ',E10.4,'  <INC)i>'  , 4X,'KA=',El0.4,15X,'b  = ',t10.4, 

, ' <1 NCH>'  ,6X» ' X3  = ' ,E10.4/) 

»AirC(6,$4)  RHO,RK 

SI  r OKMAT (//3  5X.  'CYLINDER'  , 10X, 'R= ',£12.5, '<INCH>' ,4X, ' KR  = ', El  2. 5/) 

4 Rl  IE  (6,  55) 

FOr.H  T (J5X,  'flETHCD  OF  SOLUTION  ♦MODAL*'/) 
wfi  I IE  (6,  56j  ^11 

• iAmA r (SOX,  'NORMAtlZATIOW  I Y 1 1 I = ' , E 1 0, 4///) 

«'f:t(3,57) 

*4  : (‘,5X,  'ONTA  OUTPUT  '/) 

• . . 1/6. 

•..F  A t'.oLAa  WIDTH  CF  THE  SICI 
' •»  1 4/  l2  . • RHO)  ) 
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w Jt’  INfIMXE  ShEIES 
.■i.'iAAi^=::MAX+i 
JJ  130  .1  = 1,  12 

11=.1-1 
Jt^n=i . 

*F(M.tg.1)  £2  4=2. 

2iiia1  = crf£b 

(AiS  (P(UJ*4  1-EI/2.  ) . LE.  1 .E-7)  PHIB  1 = E HI B*  1 . 00  1 
I JJ  E 1 (4)  =C0£ (41*2HI0)»  (-PI*CQS («1*FHIB1)/(( (Ml*PHIBl) ♦*2-  (PI/2.) **2 
) **2*(1./£P4) 

..wTE^tiArioii  J?  esz  (KZ)  *H1  (.'1,KZ)  *tXP  (~J*KZ*ZC)  BETWEEN  C AND  KU 
jiii.rA=  JU5HC0C  of  InE  SINOULAE  POINT  KZ=K0  IN  WHICH  THE  INTEGRAL 

.i«£  BEEN  CA..CJLATED  ANALXIICALLY 
J£LTA=0. C jO 1«K0 

J£LTA1=  ..£1  GHbOUEHOCD  OF  THE  SINGULAB  POINT  K£=KO  WHERE  THE  INTEGRAND 
i/aBIE3  a.-4PIJ..lf  ANC  'NDEbTA'  SAhPLES  HAVE  BEEN  USED. 

O£i.TA1  = 0.01^KO 
j jELTA=1 oj 

jKZb=  (DiiETrtl-OELIA)  /NDFLTA 

NSECrl=  NO.  JE  SUESECTICNS  BETKEtN  C AND  KO-DELTAl 
NSECT1=(IFI*(  (B+ZC+RHO)  ♦Ko/PI)  +2)'»‘NCYCLE 
JK21=  (KO-JELTAI ) /NSECI1 
jGECI  = N3£v:rUNDELIA+2 

j.1  = FiaST  INTEGRAL  (BETWEEN  0.  AND  KO) 

A.  1 = (O.fO.) 

JJ  20 C 1=1, NSEC  I 
aF(I. LE. N3EGT1+ 1)  GO  TO  120 
^Z=^J-D£.LTA1♦  (I -N SECT  1-2)  *EKZ2 
JKZ=DKZ2 
aJ  TJ  140 
14J  itz=  (1-1)  *JlCZl 

aF(KZ.EO.O)  rz=c. cocoi*ko 

J A Z= J K ^ 1 
l4j  waN— 1. 

aF  ( (I .DO.  1)  .JR.  (I. EC. NSECil +1) . OR.  (1 . EC  . N S ECT 1 +2)  . OR . (I. EU. NSECT) 
CIN=0.5 

i\T  = 3jRT  (KJ*KO-KZ*KZ) 

2JIE.CP=  (SIN  (KZ*E/2.  )/  (KZ*E/2.)  ) **2 *C IN » DKZ* C EXP  ( (0  . , -1  . ) * KZ *ZO) 

44AX1 =K4AX 
NJKr=RHJ*AT 

-JHPJXATaJN  JF  F4  (N)  =1./(JN  (X)  **2  + YN  (X)  »*2)  ANC  FN  (N)  = 1./ (DON  (X)  **2  ♦ 
JEh(X)*»2)  ?Jci  X = F,OKT  AND  N = 0 TO  MMAXI  ; WHERE  MBAX1  IS  A NUMBER  AFTER 
WHICH  Tiiii  :a.GTRIBUTIONS  of  FM(N)  and  FN(N)  to  THE  INFINITE  SUM 
JEC0.1E  NE.aIGIBLE.  MMAXI  IS  A FUNCIICN  OF  THE  ARGUMENT  X AND  IS  ALWAYS 
IHA«  Ja  ECUAL  TO  MHAX.  MMnXI,  FK(N)  AND  EN(N)  ARE  CALCULATED 
JY  SJBRJJTINj  FMFN (X , 4MAX1 , FM, FN) . 

0 ALL  F 4F  J (R  jKT,  K.M  AX  1 ,EM  , f N ) 

XZKTRC=  (KZ/  lKT*KO»RHO) ) 
jJ  2C0  H=1,4MAZ1 
4l=rt-1 

Hl  = (I./a1**2)  * (Fl  (4)  ♦..1*»a*K2KTRO*FN  (M)  ) 
aJvi  .1  = a1  ♦Rl  ♦2SIEXP*F  1 (M) 

Al=(2.*NO/(rA«5riO))*(I1-F1(1)*C£XP((0.,-1.)*KC*ZO)»(PI*PI/(2.*aC)) 
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j 

i 


If, 

f 


1 


L. 


j * (ol.i  J.  )/  (Kt*E/z.)  ) **2/  (2  . * (0.  577:i5fc  ony  + ALOG  (EUC*SUST 

->  I i\0/ ) ) + At.JG  (Lr,1.1A)  ) ) 

-JGPJIAIIJJi  je  12  (cETWEEi;  ZaEo  ASD  ETAKAX  ; AHERE  ETaMAX  IS  A HUMbEr 
,^Er*i2  IilE  INTEGRAND  BECCiUS  VEBX  SMALL) 

.J  = J. 

xNIrGiAriuN  IS  CAREIEL  OUT  EX  TRAPEZOILAL  RULE.  AT  FIRST  IHE  WHOLE 
a„NG2  of  IMISiaAIIOA  (C.,tIAf1AX)  IS  DtVIUEE  liJTC  TWO  SUBIMEEVALS  : 
l).,2TA1)  AGO  (tXAl ,ETAMAX)  , WHERE  ETA 1=ET AMA X/2 . . THEN  THE  NUMERICAL 
-.jilPJl'ATI  JG  OF  IHE  INTEGKAx  IS  PERFORMED  IN  THESr.  SUBINTERVALS  Wild  THE 
(J.  OF  SA.iPLES  IN  THE  FIRST  SU3INIEBVAL  TWO  TIMES  THAT  IN  THE  SECONC  ONE. 
jTAl  = 7,/  (iO”d) 

.<->BCn=  (IFIX  (S2RT  lKO*KO  + EiAl*»2)  *RaO/PI)+2)  ♦NCYCLE 
JEiAl=EXAl/NSiCn 
J iTA2  = 2.’»JiIAl 

..SEC2  2 = IFxX  ( (EIA;iaX-ETA1)/DEIA2)  +1 

NSEOr=NSEiri+NSECI2+2 

jo  joo  1=1, n sect 

xF  (I. xE. NJECT U 1)  GO  TO  220 
nrA  = SrA1+  (I-N SEC 11-2) •DEIA2 
J Ex  A=  DEI  A 2 
TO  240 

^ 2o  xIA=  (1-1)  *0EI Al 

*.F  (bPA.E'O.  3.)  EIA=0. 0001/^2 
jETA=DrIA 1 
X 40  J xN=  1 . 

„F ( (*. E^.  1)  .JR.  (I. EC. NSECT 1 4 1)  . CR.  ( I . EQ . N SEC T 1 ♦ 2)  . OH  . (I.EQ. NSECT) 

J } CIS =0.5 

PSbX=  (S.:i..i  (EIA*i;/^.  ) / (£TA*d/2.  ) ) »*24EXP  (-£TA*ZO)  *DErA»Ci:i 
Kr=S0RT(KJ*KO4ETA*»2) 

EIKIRC=  (ilA/  (RHE*KO*KI)  ) •*2 
J d A X 1 = E d « X 

CALL  FflFN  (RdO*KT,  fldAXI  ,Fd,?N) 

JJ  3(0  d=l,ddAX1 
.il  = d-1 

a1=(1/(i\T<-f.  I))*(FH(H)  -H1*M1»ETKTR0»FN  (M)  ) 
ijj  x2=IX4F  1 (1)  *i  l*tStX 
..2  = 12  *2.  ♦.<  J/  (PI4RH0) 

/l.-=(Il4().,1.)*;2)*A4B*YO/(2.*tI»PI*EHO) 

..  JS-d-xEIEAIITN  OF  IHE  PHASE  OF  Ylx 

Y..  12=Yl..*CE;<P((Ci.,1.)*(KO*SCRr  (ZC»ZC4  (R  HO*  PH  10)  **2)  ) ) 

JJdPUlAIlJ.I  OF  Td£  ACrUA..  PHASE  'PHASEY*  AND  I.OHMALIZx.1  PHASE  'I'HASNK' 
jF  YI2. 

P.i«SxY  = EAC.A.i2(UIdA0(Yl2),LREAL(YU))*1BC./ri 
PnAS.1d=DAlAA2  (DIdAG  (Yi,12)  ,DR£AL  (YN12)  ) *180./1I 

-JdPUIAIlJ;!  OF  lit  MAGNITUDE  OF  THE  Y12  IN  TERMS  OF  <«HO>  ANL  <DB>. 


^ , 3.C,  f 7.  2 , ' <0E«>'  , 10.x,  ■ DB  = ' , h12. 5 ,3X,  • KuriM.  Pii  AS  E- ' , r 7 . ^) 

j 20  P 
^.JD 

PrtOOaAW  IJ  OJ.lPBrt  IHii  MUTUAL  A Cl-ll  1 1 ANC  i:  BEIHiiEN  TWO  lEEUPICAB 
s-c:i  JA  A CYliriDER  ( UI  KOUAL  SOLUTION) 

-k~.AI.  KC,i\Z,KI,I^,KZKtro 
- 1 ) I ),Yl^,FSIr,XP,YNl2 

.(-.A-*S  ?1{4CC),  FM  (400)  ,?lj  (4CC)  ,EIi1AU,DBEAL,L)ArAlJ2 

..ti-Ur  PAEA3aIE3S: 

aJ=«XVE  iiJ.lSSS  IN  FFEE  SPACE  IS  TEHHS  OF  1/INCH 

..wYC-.E’=NU.  Jf  SUSSECIIONS  BFIASEN  ANY  TWO  SUCCESSIVE  ZEEuS  CF  IKXEGFAND 
-S  TrfAPEZJiJii^  5JLE  FOR  NUr.ERiCAL  IKIEGHAIICN 
,t»o=  SLOT  JX.iENSICN  B>A  <INCH> 

4iiO  = SACIUS  JF  CYLINDER  <INCH> 

P,ilO=  AN  JU*.AR  SEPAhAIICN  Or  I HE  SLuIS  (CENXEa  TO  CENTER)  <RADIAN> 

LJ=  SEPARAXION  CF  THE  3L0IS  IN  Z -01  EEi- 1 ION  <INCh> 

Y11=  NCSa  IE  iZaX  ICN  FACXCF, 

■11AX=  LAX^IJ.I  NO.  JF  TERMS  ntilCH  HAS  SEEN  USED  IN  CALCULAXION  CF 
iNFINIIi;  SE3Ii:3 
SMAX=36u 
21=3. 14l5^3o5 

I J=1./  (12).  -'?!) 
i\0=Z.  *EI 

^s^^  = 3. i 1 )*A0/ (2. *PI*2. 54) 

( J YCi. E=  1 C J 
4 = 0.  5 
-»  = 0.2 
tiaA  = KC'*‘A 
J a3  = N0*  3 
RH0=2  . 
a K=  KO  * R liJ 
Y 1 1 = 1 . 

7 it  I X E ( 6 , 5 J J 


J J 

r'vjRMA  X (4  C.i , 

MJIUAE  ADMITTANCE  OF  SLOTS  ON  A 

CYLINDER  V) 

*3114  (6,5)) 

JsJ 

F JnMAI  (37a, 

' AXIAL  V) 

4 RI X4  ( 6 , 5 i ) 

FREC, NO 

Format  (2  5x , 

' FREOUEKCY ' , 10X, ' F=' 

,E12. 5, ' <HZ>* 

,1CX,  ' 

K=*  , E12.  5, 

' <1/INC  .I>V) 

43aIE  (6,53) 

FJai'lni  (ICt, 

'SwOI  i-XKFHSIONS') 

4 rite  (c,53) 

.t,ANA,S,BNL 

3J 

FuuMA  r ( 1 0 i, 

•A='  ,L1C.4,'<iNCli>’  , 

4X  , ' K A=  ' , E 1 C . 

4.15X, 

•3=*  ,E10.4, 

' <I  NO  1>  ' , -tX 

, 'KB=' ,410.4/) 

« A ..  T J (t  , ; 

a HO  , R K 

F j it  4 It  1 (35a, 
* R.i  Xa  ( 6 , 5 3) 

'lYLINLEF.  ',10X,'R  = ', 

r 12. 5 , ■ <:nCH> 

',4X,  • 

KR  = ' , E 1 4 . 5/) 

> J 

F Jrt.'U  r (35t, 

'L^THCC  CF  SCDUXICN 

•MJDAi.*  V) 

« RI  (6  , 5 >) 

Y1  1 

3l> 

F jR  M«  I ( 30  A, 
*AXt£«  (c,57) 

'HORKAL^i-ATICN  |Y11| 

=' ,t1C.4///) 

j; 

F j it  .1 A . ( ’ 5 , 

Pii  10=  ?!/•♦. 

' J il  A OUT  PUT '/) 

Lt  J ^ * 
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ASJULAH  .-IDIH  CF  THE  SLOT 

{k/  (2.  *ahO)  ) 

JJhPUiAriJ.J  JF  INFINITE  StFIEb 
a.lAAl2=XRH<*1 
Jj  ICC  *1=l,rt1AX12 
n=i'i- 1 

Ir’ (.1.  1)  ;c  IC  99 

E 1 (.'i)  =COS  (Il’CHIOj  * (SIN  (Ml  ♦PHIA/2. ) / (M1*PHIA/^.)  ) **2 
j u X J ICC 
ii  E I (M)  =C  . 5 

I J\j  CJN'XlNUc 

x..IE^RAIIJN  JF  ESI(KI)  *P1  (M,KZ)  *EXP  (-J*KZ*ZC)  EETWIEN  0 AND  KO 

N£:.;rfdCUNHCOL  OF  THE  SXNOULA'n  POINT  KZ  = KO  IN  WHICH  THE  INTEGHAE 
C .»aS  dEEN  OALULATiD  ANA  L YI  IC  ALE  Y 

JELTA=1 . E-7 ♦\0 

.(SiiOr  1=NC.  IF  SAMPLES  IE  Hi  E INTERVAL  ( C NC -DELTA ) . 

NSECT1=  (IFIX  ( (S  + ZC+EHO) *KO/PI) +2)  *NCYCLE 
JXZ=(KC-Cj;LrA)/N3tCI1 
NSECr  = KS£  :ii  1 
xl  — (C*  fC*) 

-1=FIiiSI  INXEJRAL  (BETWEEN  0.  AND  KC) 

DO  ^C  0 1 = 1 , .4SHCI 
Ad=  (1-1)  *DKZ 

IF(KZ.E2.1)  KZ  = C.  C00v'11*KC 
JXE=1 . 

xF ( (I. Ev. 1)  .OS.  (I. SC. NSEC X)  ) C1N  = 0. 5 
aI  = S2S1  (K J»KC-KZ*KZ) 

xF  (AriS  (KZ»d/2. -PI/2 .)  . L£. 1 .E-H)  K Z= 1 . 0 C 000 1 * KZ 
F3IEiP=  (COS (KZ*S/2. ) / ( (KZ»3/2. ) *»2-  ( P 1/ 2 . ) • * 2) ) ** 2 *CIN * LKZ 
j*OEXP  ( (J.,-1.)  *(<Z*ZC) 

.i.UX1=E  MAX 
iOKX=  SriJ*KI 

-OMFUIAXIJN  OP  PM  (N)  =1./ (ON  (X)  ’'•2  + YN  (X)  **x)  AND  F H (N)  = 1 . / (DO N (X)  **2  + 
jY«(X)»*2)  for  X=H0KT  and  N=0  iO  MMAXI  ; WhESE  MMAXI  IS  A NUMBEF  AFTER 
WiiXC.i  Ibi  OOjTSIBUTICUS  UF  FW(N)  AND  FS(N)  TO  THE  INFINITE  SUM 
dECOME  NEGLIiXBLE.  MMAXI  IS  A FUNCTION  OF  THE  ARGUMENT  X AND  IS  ALWAYS 
xESS  THAN  0.<  ECUAL  TO  MF.AX,  MMAXI,  PR  (N)  AND  PN  (N)  ARE  CALCULAXED 
jY  SJBPGUXIN3  FRFN  (X ,MMAX1 , FM, PN) . 

CAxL  PMFN  (ROE  X, KM AX1 , FM.FN) 

JJ  ^CO  M=  1,  MMAXI 
M l = M-  1 

2Jd  x1  = Xl  ♦PN  (M)  *PSIEXP*F''  (H) 

.JMEJiATION  OP  12  (Iti«xEN  ZERO  AND  LXAMAX  ; »HERt  LTARAX  IS  A NUMBER 

.iFicA  WHICH  IHE  Integrand  beccmes  very  small) 

12  = 0. 

AiAMAX=1L.  / (20-E) 

D.ix  xNTaGaAIXON  is  CABi'Id.^  CUT  hi  TrAPELOXLAL  RULE.  Ai  FI  RSI  IHF  WOOLS 
.i.tNGE  OF  xilt-iRATION  (0.,cTAMA;<)  IS  LtVIDED  INTO  IWO  SUFINIERVALS  : 
(J.,DrA1)  A'JD  (ETAI  , SVAKAX)  , VhERt  £ lA  1 =£T  A M A X/2 . . THEN  THE  NUMERICAL 
CjMPJXA:I-N  op  IHE  Xi.TSGRAL  io  PERFOHMnD  IN  THESE  SUEINTERVALS  WITH  THE 
..0,  jF  JA-lFLPj  IN  Tl.E  FIRST- SUUINXSPVAL  InC  TIMES  XilAT  IN  THE  SFCON’D  ONE. 
cXrt1  = 7./  CJ-Jj 

. jwC  X 1=  ( I FI  X (S2R  T (.^C  ♦ ETA  1 ‘‘x)  » RiiO/PI ) *2)  *NC  YCLi. 

Jwl/t1  - 1 1'  f 1/  iS.CII 
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A2  = 2. ♦ J£IA  1 

Ia  { {ETAMA»£rA1)/CEIA2)  ♦I 
-.iiiCA=N3t  JI  UNSECI2  + 2 
ju  it  0 1=1,  NScCI 
xi?  (I.  LE.  Ei3C;  1+1)  GC  10  220 
^EA  = aTA1+ (I-iiECrl-2) *DETA2 
JiTA=D£IA2 
iJ  IJ  2<4C 

^2j  22A=  (1-1)  *0ErAl 

I x£(EaA«E^»0*)  £TA~0«C0C1/A 

i,  J2rA=C£IAl 

.)  24iv^l1~1« 

‘■-I  -r'(  (^.E2.  1)  .02.  (I.EC.NSECl'1+1)  .CE.  ( 1 . 2C  . NS  ECU  *2)  .03.  (I.EO.NSECT) 

'*!  ^)  1-23=3.5 

2iEX= (CCSi  (2rA*2/2. ) / ( (EIA*B/2. ) **2+  (El/2. ) **2) ) **2*DEIA*CIN 
^•'EXP  (-EIA*20) 
a2=S23T  (K0*KJ  + EIA**2) 

.liAX1=KMAX 

JAxL  FMFN  (2iO*KI,  flMAXl,FM,FN) 

jj  3C0  a=i,.'i.uxi 

JJj  x2  = I2+Fi  (.1)  ♦PiEX*f1  (M) 

2 12=  (I1»  ( ).  , 1.  ) *12)  *A*B*YO/  (P2*KC*P.H0**2) 

C J03HAL22A22C3  0?  IHE  PHASE  OF  Y12 

YN12=Y1  2^CS:tE  ( (C.  ,1 . ) • (KC’SCET  (ZC*ZC+  ( S HO* PHIO ) **2)  ) ) 

C - jaEJIAIIGN  JF  IHE  ACTUAL  PHASE  'PHASEY'  AND  NOEMALIZEt  PHASE  'F.HASLK* 

C jE  YIa. 

Ei1ASEY=DArAN2  (UIHAG  (Y12)  ,LREAL  (Y12)  ) *180./PZ 
EiiHSSn=  jAi'A.12  (DiHAG  (YN12)  ,DSEAL  ( Y N 1 2}  ) * 1 SO  . / PI 
C -JlIPJXAriON  JF  IHE  BAGNITUEE  OF  IHE  Y12  IK  lEEMS  OF  <KHO>  AND  <DB>. 

...iPY=CLAEi  ( Y 1 2) 

.i.lPYOa=ALJi  n (AMPY/A3S  (Y11)  ) *20. 

.tPr.IK  = KC*  ,iHJ*£HIO 
-ON=KO*ZC 

P.ixOD  = Prll  J*  1i  1./PI 

2.X12E  (b,4Lj)  PHIOL,EPriIK,ZO,ZOK,AMPY, PHASEY, AMPYDB,PHASNN 
Jj  FOEHAT  (2CX//1  IX,  • PH10='  ,F7. 2,  '<LEG>  * ,4X  , • K*E*PHIO=  • , El  0.  4 

, IbX,  *Z0='  , El  1. 4,  *<INCH>'  ,4X,'K*Z0='  ,i1C.4/10X,  • Y12=*  , E12.5, 

^ ' <M  H0>'  , it,  ?7. 2 ,'  <DEG>'  , 10X,  'EB=»  ,F.12.  5, 3X,  ‘NORM.  PH  AS  E= ' , F7 . 2) 
ii’OP 
-HD 

iJBEJUTLNE  F4FN  (X,N,FM,FN) 

2EAL*i  CU.II  (-*CC)  , CUM2  (400) 

2-AL*8  FJ  (401)  ,XB,bSSY  (400)  , FM  (400)  ,FN  (400) 

PI  = i. 14  1 5 12  ji 

A i = X 

xF(X-0. 1)  1 10,20 
10  jai’ixJG  = AiJG  (.(/a;.)  *0. 5772156(349 
A 2=  X*  X 
A J = \2  *X 

rA4=X*Xi 

A J=X«X4 

j JJ  Y I = A.  * ( 1A  UOG*  (1  . - X2/4.  ♦y.4/6  4. ) ♦X2/4.  -3.  *XC/12S. ) /P  I 
J3j  Y2=-2.  / (PI  *X)  +2.  * (GAllL^G*  (X/2.  -X3/16.  ♦X5/364. ) -X/4.  ♦ 1 . 25  * X3/ 1 6 . 
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<i  ( 


^-3,3j33J*X5/7o6.)/PI 
3 j r j 2z 

2j  -«i,L  BEJY  (X,  J,BS3Y1,IHP) 

3ALL  ilESY  (X,  1 ,&SSY2,IhK) 

33  JOnriNL’E 

3S3Y  11) =E3ii 1 
J3SY (2) =ES3Y2 
JJS3Y 1=-E33Y (2) 

2=1 

^=1^1 

J33Y  (I+  1)  =2. ♦ (1-1) •ESSY  (I) /X-BS3Y  (1-1) 

*>33YI  1 = aS3Y  C»1) 

2F  (Ad3  (BS3Y2  1)  . GE.1 . ElO)  3C  TC  100 
3J  TO  63 
1 jJ  JMAX=I+1 

I?  (M3AX.  GE.  i)  BEAX=.M 
3l=H-1 

GAEE  aSLJE  (Xd,PJ,  tit1AX  + 1 ,0.C00,7,IERR,DU(n,  CUH2) 

JFJ1=-FJ  (2) 

Fii  ( 1)  =1.  / (BS3Y  ( 1)  ••2  + f J (1)  ♦♦2) 

FtJ(1)=1./(D333Yl**2+DFJl*»2) 

JJ  2C0  I=),N1 
^F (I. 3E. K^AX)  GC  TO  250 
JdS3Y=BSSY  (I)  -I*BSS Y (I*1) /X 
JFJ=FJ  (i)  -I*FJ  (!♦  1)  /XE 
Fil(i*1)=1./(3  33Y  (1+1)  **2  + FJ  (1  + 1)  **2) 

F.J(I*1)  =1./  (0BS:y**2  + DFJ**2) 

2Jv(  JdUIIdUE 
2iJ  JJBTISUi 
i=NHAX 
ddl'JXN 
3dd 

G dJSSOUTIEi  'diSY' 

C PUBfOSd 

C COHPJIE  IHE  Y BESSEL  FDKCTION  FOE  A GIVEN  ABGOHENT  ANE  ORDER 

G USAGE 

CSLu  3E3Y  (X,N,BY,IES) 

oEscaipnod  of  parameters 

X -THH  ARGUMENT  CF  THE  Y BESSEL  FUNCTION  DLSIREt 
N -TdE  ORDER  OP  THE  V BESSEL  FUNCTICN  DESIRED 
BY  -THE  RESULTANT  Y BESSEL  FURCTICN 
lER -RESULTANT  ERROR  CODE  NhERE 
223=0  NC  ERROR 
IE3=1  N IS  NEGATIVE 
2S3=2  X IS  NEGATIVE  OB  ZERO 
IE3=3  BY  HAS  EXCEEDED  MAGNITUDE  CF  10**70 

REMARKS 

VEFY  SJaLL  VALUES  Or  X MAY  GAUSt  THE  RANGE  OF  THE  LIBRARY 
fUKGriON  ALCG  TO  PE  EXCEEDED 
>■  MUSI  BE  GBEAIcR  IHAN  ZERO 
N KJ3I  dE  GRuAIER  THAN  03  EOUAI  TO  ZERO 


n o cj  ( i ( . ( . r.  « j cj  c/  < - r,  f j n f f . r j fi 


oJBROUriNES  ANC  FUNCTION  SUBFGOGBANS  BEgUIfifC 
NONE 


ilirHCD 

bcCJRBENCi:  HEIAIION  AND  PCiYllOMIAL  AFPGOXIHATION  TECHNIQUE 
AS  CSSC3IEEC  BY  A . J . H . HITCHCOC K , ' POL YNOHI AL  APPRCXIMATIONS 
TO  SE3SEL  FUNCTIONS  OF  OBOEB  Z£HO  AND  ONE  AND  TO  BELATED 
F JNOIIO.'JS  ' , M.T.A.C.,  V.  1 1 , 1 95  "J,  P P.  66-8b  , AND  G.  N.  HATSON, 
■A  IHEAIISS  ON  THE  IHEOBY  OF  BESSEL  FUNCTIONS',  CAHBBICGE 
USIVEBSIIY  PRESS,  1950,  P.  62 


jUBROUTINe  E E3 Y (X ,N , B Y , I ER) 

CHECK  FOE  2SHJBS  IN  N AND  X 
IF(N)  180,10,10 

1 u i.  2B-=0 

.F (i)  190,19  0,20 

jBANCH  if  X LESS  IHAN  OR  ECUAL  U 

2j  iF  (X-4. C) 40,4 J, 30 

COMPUTE  YO  ANI  Y1  FOR  X GREATER  THAN  4 

Jj  11=4. 0/X 
C2=ri*T1 

F0=  ( ( ( (-.0000 03 70  43 *T24. 0000 173 565) *T2-. 0000487613) *T2 
i 4.00017343) *Ti-.C017530b2) *T24.3989423 
J0= (( ((. COO  00 32312*12-. 0000142078) *124.0000342468) *T2 
t -.C0C38o97i1) *124.0004564324) *T2- . 01246694 
P1= ( ( ((. 000 0042 41 4*T2-. 0000 200920) *124.0000560759) *T2 
j -.000^23203) *12 4. 002921826) *T24. 39894^3 
Ql=  ((  ( (-.  OC  0 0 03 1>5  94  *T2  4. 00  00  162  2)  *T2-.  0 0003  9 6708)  *T2 
t 4.U0C1C J4?4l) *12-. 0006390400) *124 . 037 40064 
<i  = 2.  O/SCRI  (X) 
j=A*n 

C=X-. 7853382 

X0  = A*PC*SIJ  (C)  4 E*  to*  COS  (C) 

Y1=-A*?1*C0S(C) 4B*Q1*3IN(C) 
iJ  TO  90 

COSPUIE  YC  AND  Y1  FCK  X LESS  THAN  CB  EQUAL  TO  4 

•♦0 

a4=Xa*>;.( 

C = ALJG  (aX)  4. 577a1 57 
s Jrt  = 0 . 

2 8a«=i 
M=T 

«.JJ  7j  L=1,1j 
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n Ci 


(L-1) 5C,60,50 
50  jdM=jUH*1./FI,OAI  (1-1) 
oJ  gL=L 

i'5  = I-3Ufl 

rEhM=  (TEE. 1*  (-X2)/FL**2)  • (1,-1./  (FL*TS)  ) 

7j  X0=Xu+T£B.1 

2EHtl  = XX*(I-.5) 

j0M=0. 

f 1=TEEM 

JO  80  1=2,16 

j JM  = 5 UK* 1. /FLOAT  (1-1) 

Fi.=  l 

FL1=F1-1. 

i3=I-SQH 

2EEH=  (lEE  1* (-X2)/  (FlI^Fl) ) * ( (TS-.5/F1)/  (TS  + . E/Fll) ) 
Jj  X1  = n+TcE.1 
^ i^=,  63o6193 
Xj=PI2*X0 
X 1=-PI2/X*PI2*Y1 

JdECX  iF  JJIX  YO  CE  T1  IS  CESIRFE 

iJ  ^F(JJ-I)  1C3,  1DC,  130 

iETURS  £i:8£3  YO  CR  Yl  AS  EEQDIfiEO 

1 Jj  IF  (S)  110,  120,  lie 
1 U JY=Y1 

JJ  TO  17C 
Izj  JY=Y0 

JO  TO  17C 

P::.EFJE.i  ££OUBa£NCB  OPERATIONS  TO  FINC  YN  (X) 

13j  YA=Y0 
/J=Y1 
X=1 

14J  T=F10AT (2*K) /X 
YO=T*YB-YA 

IF (A JS  (YC) - 1. CE70)  145, 145, 141 
l4l  IEH=3 
BET JB  N 
145  E=K4l 

iF (K-N)  150,  160,  150 
13J  Xii=Y3 
X8=YC 
JJ  IJ  14C 
loO  ilY=YC 
17j  BEUaj 
1d««  iBB=1 
BBTUBN 
1 iu  i.dR=2 
BSTUaN 
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SOBROUTIN’'  BSLJZfr  , FJ  , NBAX  , X , ND  , lERR  , FJUPEX  , RP) 


THIS  IS  ONE  OF  THREE  ROOTIHES.  "BSLJZ",  "BSLIZ",  AND  "BSCJZ", 
BASED  ON  ALSORTTHH  236  FHOH  "COBHD N IC ATI ONS  OF  THE  A.C.B.", 
ADGUST  196U.  -^RTS  ONE  EVALUATES  THE  BESSEL  FONCTIONS  OF  THE 
FIRST  KIND  FOR  REAL  ORDERS  AND  NON-NEGATIFE  REAL  ARGUMENTS. 


THE  PAHAHETFRS  ARE  DESCRIBED  AS  FOLLOHS,  WITH  AND 

"(I/O)"  INDICATING,  RESPECTIFELT,  THAT  A PARAHETER  IS  TO  BE  SET  ON 
ENTRY,  WILL  BE  SE"’  BY  THE  NOUTINB,  OB  BOTH  : 


ALL  PARAMETERS  ETCEPT  "ND"  , "lERH"  , "NMAX"  ARE 

••*  DOOBLE  PPSCTSION. 

(I)  X — THE  (NON-NEGATIVE)  ARGUMENT  TO  THE  BESSEL  FONCTIONS. 

(Oj  PJ  --  AN  ARRAY  IN  WHICH  THE  VALUES  OF  THE  BESSEL  FONCTIONS 

ARE  STORED,  AS  FOLLOWS:  LET  J(X;B)  DENOTE  THE  VALUE 
OP  THE  BESSEL  FUNCTION  OF  ORDER  B WITH  ARGUMENT  X. 
THEN,  FOR  I = 1 TO  ABS(NMAX)+1, 

FJ(I)  = J(X;A  ♦ (1-1)  •SIGN  (NMAX)  ) , 

(I)  NMAX  --  REFER  TO  "FJ". 

Id  a --  REFER  TO  "FJ".  NORMALLY,  0 <=  A < 1,  BUT  THE  ALGOR- 

ITHM WORKS,  WITH  SOME  LOSS  OF  ACCURACY,  FOR  A >=  1 . 
SEE  THE  PROGRAM  NOTES  BELOW. 

(I)  ND  --  THIS  GIVES  THE  NUMBER  OF  SIGNIFICANT  FIGURES  OF 

ACCURACY  DESIRED  IN  THE  FUNCTION  VALUES. 

(0)  lERR  — THIS  IS  AN  ERROR  FLAG  WHICH  IS  SET  TO  0 IP  THE 

INPUT  PARAMETERS  ARE  OKAY,  AND  TO  SOME  POSITIVE 
VALUE  IF  ONE  OF  THE  PARAMETERS  IS  INVALID.  REFER 
TO  THE  ERROR  EXITS  AT  THE  END  OF  THE  CODE  FOR  A 
DETAILED  LIST  OF  THE  VALUES  OF  lERR. 

(0)  FJAPRX  --  A SCRATCH  ARRAY  USED  BY  THE  ROUTINE.  IT  MUST  HAVE 

AT  LEAST  ABSjNHAXWI  ENTRIES. 

(0)  RR  --  ANOTHER  SCRATCH  ARRAY.  IT  TOO  MUST  HAVE  AT  LEAST 

ABS (NMAX) +1  ENTRIES 

OTHER  ROUTINES  CALLED:  ( • INDICATES  A LOCAL  ROUTINE  ) 

• NBSOn  --  INVERSE  FUNCTION  OF  X*LOG  (X) 

» UNDBRZ  — ROUTINE  TO  CONTROL  UNDERFLOW  INTERRUPTS. 

D3AHMA  — DOOBLE  PRECISION  GAMMA  FUNCTION. 

DLOG  --  DOUBLE  PRECISION  LOGARITHM 

DABS  — ABSOLUTE  VALUE 

MOD  — REMAINDER 

DMAX1  --  MAXIMUM  OP  2 REALS 

NOTES: 

THE  METHOD  OP  COMPUTATION  IS  A VARIANT  OP  THE  BACKWARD 
RECURRENCE  ALGORITHM  OF  J.C.P. MILLER  (REFERENCE  1).  THE 
PURPORTED  ACCURACY  IS  OBTAINED  BY  A JUDICIOUS  SELECTION 
OF  THE  INITIAL  VALUE  "ND"  OF  THE  RECURSION  INDEX  (REP- 
RESEN'"ED  IN  THE  CODE  BY  THE  VARIABLE  "XNU"),  TOGETHER 
WITH  AT  LEAST  ONE  REPETITION  OF  THE  RECURSION  WITH  "NU" 

REPLACED  BY  "N0"*5.  NEAR  A ZERO  OF  ONE  OF  THE  BESSEL 
PDNCTTONS,  THE  ACCURACY  OF  THAT  PARTICULAR  BESSEL  FUNCTION 
MAY  DETERIORATE  TO  LESS  THAN  "ND"  SIGNIFICANT  DIGITS.  THE 
ALGORITHM  IS  MOST  EFFICIENT  WHEN  X IS  SHALL  OR  MODERATELY 
LARGE. 


THE  ABOVE  PARAGRAPH  IS  TAKEN  FROM  GADTSCHI'S  PRESENTATION 
OF  ALGORI’^HM  236  IN  C.A.C.M.  THE  SELECTION  OF  THE  INITIAL 
"NO"  IS  DONE  WITH  THE  AID  OF  THE  FUNCTION  N3S01Z,  ALSO 
BY  GAUTSCHI  (AND  CALLED  "T"  BY  HIM).  IN  THIS  CODE,  THE 
FOLLOWING  SP’^CIAL  CASES  HAVE  BEEN  ADDED: 

A,  r=0  WHEN  NMAX  > D OR  A=0 

B,  A=0  AND  NMAX  < 0 

C,  A >=  1 : THE  ALGORITHM  WORKS  IN  THIS  CASE,  BUT  THE 

INITIAL  CHOICE  OF  "NU"  IS  NO  LONGER 
OPTIMAL,  AND  SOME  ACCURACY  IS  LOST,  SIMPLE 
TESTS  INDICATE  THAT  ONLY  A FEW  DECIMAL 
PLACES  ARE  SACRIFICED  AT  HORST,  A LIMIT  OF 
"ABIG"  IS  PLACED  ON  A TO  AVOID  OVERFLOW  IN 
THF  GAMMA  FUNCTION.  TO  AVOID  COMPLICATIONS, 
NMAX  IS  REQUIRED  TO  BE  NON-NEGATIVE  IF  A > 
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REFEEEHCHS: 

1.  rjS'lT'^CHI,  W.  "'>Ecn'’SI  VE  COMPOTATION  DF  SPECIAL  FUNCTIONS", 
UNTV''PSirY  OP  1ICHIGAN  ENGINEERING  SUMMER  CONFEP- 
EVCES,  VUrEPTCAL  ANALYSIS,  1963. 


PR  (1) 


* 


*****««»*• 

SUBROOCINE  PSlJZ(r  , F.l  , NMAX 
IMPLICIT  ?-"AL*q  (A-H-0-Z1 
DOUBLE  PPSCTSrOV  NBS01Z 
DIMENSION  FJ(1)  , FJAPRX(I) 

LOGICAL  NEVFN  , AFLAC 

DATA  0!IE/1d5/  , 

TPN/1OD0/  , 

CP/-*.  3P91DC/  , 

ZEPO/^D"/  , ABic/55D0/ 

ALEPH /ZUFFFFFFFFFFFFFFF/ 
C'^/ZUFAUFOOUOOOOP'^UO/ 


ND 


TWO/2DO/ 

3MALL/1D-15/ 

C3/2.3026D0/ 


lERP  , FJAPRX  , PP) 


HALF/. 5D0/ 
CV.  73576D0/ 
C4/1. 3863DO/ 
THOP5/2. 5D0/ 
FOOH/4DD/ 


k • ** 


INITIALIZE  THE  '•p'>OG  PARAMETER  , TORN  UNDERFLOW  OFF  , AND  CHECK 
THE  PARAMSTIRS  ■'OP  VALIDITY  AND  FOR  THESE  SPECIAL  CASES: 

A.  X'P  WITH  NMAX  > 0 OR  A=0 

B.  A = ''  AND  NMAX  < 0 


THE  COD^  DEL:BSRA~’='LY  avoids  TESTING  MORE  THEN  ONE  THING  IN  EACH 
LOGICAL  "IF"  B'lLOW  BECAUSE  OF  I.B.M.  FORTRAN  INEFFICIENCY  IN  THIS 
PEGAPD. 


IF  A>1,  NMA’'  MUST  NOT  EE  NEGATIVE. 

««***« 


IER°  = '' 

CALL  UNDF'’Z  (’OF'''  , 


!'• 

A 

.LT. 

ZERO  1 

IF 

A 

. GT. 

A BIG' 

I F 

V 

. LT. 

ZE"0'i 

IF 

NMA'' 

. GE. 

I’' 

A 

. EO. 

ZERO) 

IF 

A 

. LE. 

SMALL 

A 

. GF. 

ONE 

IF 

X 

. GT. 

Z'’SO) 

I'’ 

NMAX 

. GF. 

IF 

. GT. 

■^FRO) 

****«:« 


SAVE) 

GOTO  999 
GOTO  998 
GOTO  997 
GOTO  10 
GOTO  10 
) GOTO  996 
) GOTO  99U 
GOTO  ur 
GOTO  20 
GOTO  995 


IF  UMAX  < 0,  NMAXT  IS  SET  HERE  SO  THAT  ONLY  J(X:A)  IS  CALCULATED. 
THE  LOOP  FOLLOWING  STMEMENI  8C0  THEN  CALCULATES  THE  REMAINING 
FUNCTIONS  BY  A SIMPLE  PECUPRENCF. 

I"  A=0,  NMAXT  IS  SET  SO  THAT  J(X;A+N),  N=0 NMAX  ARE 
CALCULATED;  THE  CODE  AFTER  8U0  THEN  REVERSES  THE  SIGN  OF  EVERT 
OTHER  ONE. 


i 


f 


W''  FIRST  HANDLE  THE  CASE 


2C 


N'EMP  = TABS  (NMAX)  ♦ 
DO  30  I = 1.N'"EMP 
FJ(I)  = ZF®0 


x=o. 

1 


I \ 


1 T 


(••««*** 

4''  A='t.A’.  = (A  .'0.  .AND.  (Nr*AX  . LT. 

N'iA:r  = viAy 

.lt.  'i-ArT  - i 

fj  T r>  = * s Y * ( >;  •<  » Y ♦ 1 1 ) 
lOTO  f"' 

N *1  A ^ - '.’"AY 

N r F Y P - ‘i  “ ^ V ~ ♦ ’ 

Sr  E^sijN  - -FV*  ' [-'.’D) /2 

DO  ■lO  : = 1, 

9F  FJAr97(T)  = " 

S'i<i'=  Y/-vi■^)  ••  A1y:\  (0»,f+^) 

D1  = ♦ -;i 


^ = -T  •:>  O 

r='(.)MAy'  .0'.  0)  » = S-YA'--  • NBSCIZ  (HALP*D1/N:1AXT) 

S = C2  • Y . N3S01Z (Cl *D1/X) 

(;*•**•« 

c THE  ?f:jfs:''s  r'n'^x  "vh"  :s  deliberately  calculated  as  a floating 

C POINT  NUY.pYn  ?»-:i7P  ■"  «f.  ^ ah  INTEGEP,  AND  ALL  COMPARISONS  WITH  IT 
C ARE  DONE  AS  FLOATING  PCINT  COMPARISONS. 

^ * * * * « * 

v>rj  = -AH?  + DM  AT  1(0,3) 

XL  I MIT  o XNU/2 
r«OA  = A * A 
XN  = ■”’0 

FL  o 

3«  « **«  * 

C TffF  OUTER  :'"5A’T^,V  LOOP  STARTS  HERE. 

r 


C 

C THE  FOLLO'rf*'-r;  tt^d  T3  nONE  ENTIRELY  IN  FLOATING  POINT  FOR 
C EF'-ICIENCY. 

(3*  • **  * * 

200  XN  = XH  + ^N" 

FL  = FL  * (XN  ♦ A)/(XN  + ONE) 

I'='{Y.N  .LT.  2C0 

OLDFL  = FL 
OLOXS  = XN 
C 

N = 2 * Y «; 

XN  = N 

NSVFN  = .'on-. 

R = 7,  ERO 

s = z X^O 

"EXPI  '=  ’'W'^/Y 

^•* • * * ♦ 

C IN  "HY  POLLONTNG  LOOP-  Ti{  F SUCCESSIVE  VALUES  OF  "R"  ARE  PARTIAL 
C FRACTIONS  0®  A CONTINUED  FRACTION, 

(;*•**** 

3C0  DFNOM  = ?E"P1  • (>  + XN)  - R 

IF  (0  ADR  (D'"NO’‘)  .LF.  SMALL)  DFNOS  = DENOM  ♦ SMALL 
P = CNF/DRN''M 


(1|  = SUB/(OH*:  ♦ S) 
(NMAXT  .EO.  Oi  GOTO 


PLH3DA  = ZERO 

IF  (.NOT,  NEVER)  GOTO  400 

FL  = FL  * (XN  ♦ T«JO)/(XN  ♦ TWOA) 

FLB3DA  = FL  * (XN  ♦ i) 

433  S = R * (FLBBDA  ♦ S) 

I? (N  . LE.  NBAXT)  RP (N)  = R 
N = N - 1 
XN  = XN  - ONE 
NEVER  = .NOT.  NEVEH 
IP(N  . GE.  1)  GOTO  300 
0««*«** 

FJ  (1|  = SOB/(ONE  ♦ S) 

IF  (NBAXT  .EO.  0)  GOTO  603 
DO  503  N = 1, NBAXT 
500  PJ(N4^1)  = RR(N)  * FJ(N) 

C THE  LATEST  A PPROTTB ATIONS  ARE  CHECKED  FOR  IHPROVEHENT: 

600  DO  900  N = 1,NTEHP 

IP (DABS (FJ (N)  - FJAPRX(N))  . LE.  DABS  (P J ( N)  ) *EPSLON)  GOTO  800 

DO  700  rt  = 1. NTEBP 
700  FJAPRT  (B)  = FJ  (B) 

XN  = OLDXN 

FL  = OLDFL 

XLIBIT  = XLIBIT  ♦ THOP5 
GOTO  200 
800  CONTINDF 

I?(NBAX  .GE,  0)  GOTO  1003 
C****** 

C IF  NBAX<0,  WE  HAVE  FINISHED  OBTAINING  J(X;A(  , AND  NOW 
C ITERATE  TO  FIND  ALL  THE  DESIRED  FONCTIONS. 

C 

C FIRST  WE  CHECK  FOR  THE  SPECIAL  CASE  A=0, 

IF(.NOT.  AFLAG)  GOTO  850 
NBAXT  = -NHAX  ♦ 1 
DO  820  N = 2, NBAXT, 2 
820  FJ  (N)  = - FJ  (N) 

goto  1000 

C****** 

950  FJ(2)  = TWO  * A * FJJII^a  - FJ(2) 

IF  (NHAX  .EQ,  -1)  GOTO  1000 

C THE  FOLLOWING  CODE  IS  A RENDITION  OF  THE  LOOP 
C DO  900  N = 2, NBAXT 

C 900  FJ(N*1)  = (2/X)  ♦ (A-N)  *FJ  (N)  - FJ(N-I) 

C 

C WITH  OVERFLOW  DETECTION.  AS  SOON  AS  THE  NOHBERS  GET  TOO  BIG,  THEY 
C ARB  SCALED  DOWN  (BY  A POWER  OF  THE  HACHINE  BASE,  SO  AS  TO  AVOID 
C LOSS  OF  PRECISION)  AND  THE  CALCULATION  CONTINUES.  A SEPARATE  LOOP 
C TRANSFORHS  THE  SCALED  VALDES  TO  THE  CORRECT  OUTPUT  VALUES,  SETTING 
C TOO-L.\HGE  ONES  TO  PLUS  OR  HINDS  INFINITY. 

C****** 

NBAXT  = -NHAX  ♦ 1 


AND  NOW 


FJNMI  = FJ(25 
THBP1  = TWO/X 
OVER  = ZERO 
XNBi  = TWO 
C 

DO  830  B = 3.NBXXT 

FJN  = TEBPI  * (A  - XHH1)  • FJHHl  - FJNB2 

FJ(f32  = FJNB1 

FJHB1  = FJH 

FJ(B)  = FJN 

XNBi  = XNB1  ♦ ONE 

RR (Nl  = OVER 

IP/OX0S(FJ»)  ,LT.  C5)  GOTO  880 
OVER  = OVER  ♦ ONE 
FJNBI  = FJNB1/C5 
FJHB2  = FJNB2/C5 
880  CONTINUE 
C 

lEfNBAXT  .LE.  3)  GOTO  1000 
OVER  = ZERO 
SCALE  = ONE 

DO  900  N = U.NBATT 

IF (OVER  ,LT.  FOUR)  GOTO  890 

FJ(N>  = DSIGN (ALEPH,FJ (N)) 

GOTO  900 

890  IF(RR(N)  ,GT.  OVERI  SCALE  * SCALE  • C5 
FJlN)  = FJ(N)  * scale 
OVER  = RR(N) 

900  CONTINUE 
GOTO  1000 
C****** 

C ERROR  EXITS  FOLLOW,  HEARINGS  OF  THE  EXIT  VALDES  OF  "lERR"  ARE 
C 0 ; NO  ERROR 

C 1 : A < 0 

C 2 ; A > ABIG 

C 3 : X < 0 

C U : NBAX  < 0 AND  0 < A < SHALL 

C 5 : T=0,  NHAX  < 0.  AND  A > 0 

C 6 : NBAX  < 0 AND  A >*  1 
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S 

lERR 
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1 

995 

lERR 

= 

lERR 

♦ 

1 
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lERR 
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1 
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lERR 
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999 
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3 

IFRR 

4 

1 

1000 

CALL 

UNDERZ ( 

1 

S' 

RETURN 

END 

